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ABSTRACT 


TkU xzpoHt oi tfuiee. poAt^i 

VoJit I * Stmdy Statz Vzcoupling 
?(VU II : Stability and VzAlgn 

Pant III; Application to STOL AiAcaaf^t 

Pant I paojizntii a con&t^ctivz cAitzxion decoupling 
the -Steady states o^ a limaa time-inooAiant multivoJiiablz 4^4- 
tm. Thu QAltzxion con&ut6 o<j a 4et o^ inequalities Mhick, 
(A)kcn satisfied, tAiiiZ cojuse the 4itzady states oj5 a 6y4>tm to be 
decoupled. It tuJim out that puae integAotou in the loop6 
play an important aole. Stability analysis and a new design 
technique ^oa 4uck 6 y^ terns cuie given in Paat II. A new and 

I * 

simple connection between single-loop and multivaAiable cases 
is iound. This makes possible the application o^ the existing 
Single-loop methods to multivariable cases. These results are 
then applied in Part III to the compensation design ior NASA 
STOL C-SA aircraft. Both steady-state decoupling and stability 
are justified through computer similations. 
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NOMENCLATURH 


1 . 

r 

nxl input vector 

2, 

y 

nxl output vector 

3. 

H(s) 

closed-loop transfer function matrix 

4, 

Gp(s) 

nxm plant matrix 

5. 

G,(s) 

mxn compensator matrix 

6. 

T 

P 

nxm type number matrix of the nxm plant Gp(s) 

7, 

T 

mxn type number matrix of the mxn compensator 


c 

Gc(s) 

8. 

P 

poles of any transfer function g(s) 


g 

9. 


zeros of any transfer function g(s) 

10. 

Max{ * * ’ } 

maximum value among all the elements in the 
brachets 

11. 

LCD{ • * * } 

least common denominator of the elements in 
the brackets 

12. 

LCM{‘ • ‘ } 

least common multiplier of the elements in 
the brackets 

13. 

Vii 

1 

cofactor of the ijth element of 

14. 

(I*GpGc)i^ 

cofactor of the ijth element of I+GpG^ 

15. 

detd + GpGj.) 

determinant of the matrix I+GpG^ 

16. 


minor of matrix G formed from rows 



and columns j j ‘ , j 


m 
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1, INTRODUCTION 

Considerable research has been done on the decoupling of linear multi- 
variable systems (e.g. see (1)). Such decoupling, referred to as dz~ 

COuptAJig in this report, requires the system to be characterized by a non- 
singular, diagonal transfer function matrix, and in general, linear state 
variable feedbacks have been employed. 

The advantage of total decoupling is obvious, however, due to the re- 
striction of having a diagonal transfer function matrix, less freedom should 
be expected when stability of the system is concerned. 

This loss of freedom can be recovered to some extent by requiring only 
the steady states to be decoupled. Loosely speaking, a steady-state decoupled 
system is one in which changes in each input (i-th) are reflected in a cor- 
responding output and only that output, when steady state is reached. Thus, 
different from total decoupling described above, mutual interactions are 
allowed during the transient period (but only during this period). 

Necessary and sufficient conditions for decoupling the steady states of 
a system via linear state variable feedback were obtained by Wolovich (2) . 

His result, in terms of transfer function matrix representation is as 


follows: 
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A system characterized by an (n x m) proper rational. 

transfer function matrix, G (s) having no poles at the 

P 

origin (s = 0) , can be steady- state decoupled (via 

linear state variable feedback or perhaps some other 
less ambitious scheme) if and only if 

p(G (0)) = n (1-1) 

where p(G (0)) denotes rank of the matrix G (s) as 
p P 

s approaches zero . 

However, it is found that if classical cascade feedback compensation 
other than linear state variable feedback is used, the rank condition (1.1) 
is no more necessary. Furthermore, the precluded poles at the origin are 
allowed. Actually, such poles are very helpful for decoupling the steady 
states of a multivariable system. Therefore, significant advantages over 
the linear state variable approach can be obtained through classical feed- 
back configuration which then is obviously not "less ambitious". 

The constructive criterion for steady-state decoupling will be derived 
in this part of the report. It will be shown that this criterion consists 
of n(n-l) inequalities ( n is the number of outputs of the given plant), 
with the type numbers of the compensator transfer functions as unknowns. 

These unknowns are chosen to satisfy the inequalities and hence achieve a 
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steady-state decoupling scheme. Fundamental mathematical relations are de- 
rived in Chapter 2. Two simple applications for 2x2 and 3x3 cases 
are given in Chapter 3. Finally, the general case is considered in Chapter 
4. Direct comparison of the result to that of the state variable approach 
is included in Chapter 9, which marks the end of this report. 

The research reported herein was included in the Jen-Yen Huang M.S.E.E, 
thesis at the Department of Electrical Engineering and Computer Science, 
University of Stanta Clara, Santa Clara, California. The thesis was super- 
vised by G. J. Thaler, U.S. Naval Postgraduate School, Monterey, California. 
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2. FUNDAMENTAL RESULTS 

The system under consideration in this thesis is shown below in 
Figure 2.1: 


^ y 


FIGURE 2.1 

Where Gp(nxm) characterizes the given m-input, n-output plant, 
Gj*(inxn) is the n-input, m-output compensator to be designed. 

N unity feedbacks are used and complete controllability and 
observability are assumed (3), (4), to assure the complete 
description of the system by transfer function matrices. r, y 
are the nxl input and output vectors, respectively. 

Let H(s) * ^^ij^®^^nxn closed loop proper transfer function 

matrix, then by the above assumption, it characterizes the system 
completely, and we have: 




COMPENSATOR 

N. 

PLANT 


; 


G^Cmxn) 


G (nxm) 
V 



5 


y(s) = HCs) t(s) 


or 


n 

(s)= ^ (s) (s) 

j = l ^ ^ 

n 

= h..(s)r.(s) + I hi.(s)r.(s) (2.1) 

11 1 XJ J 


By (2.1) and the Final Value Theorem, we have: 




lira y. (t) = lim sy^(s) 
t -^°° ^ $"►0 


= lim s h^.(s)*r^(s) 

S -»“0 


+ lim 

S**’0 


n 

S I h..(s)T.(s) (2.2) 
J=1 ^ J 


i»l, ,n 


Then the following formal definition can be given: 


DEFINITION: 

A system with the transfer function matrix H(s) is steady-state 
decoupled if and only if it is asymptotically stable^ and 


lim 

s-^o 


n 

Z s h- . (s) T . (s) 
j=i ^ 


O 


for all i=! 1 , * • * , n 


(2.3) 


^i.e. , all the poles of the closed loop system lie in the open 
left half plane (Re(s)<o). This guarantees the application of the 
final value theorem. 
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For systems as shown in Figure 2,1, it is well known that the 
closed loop transfer function roa tr ix H ^ can be expressed as 


H = (I + G G G G 
pc pc 


(2.4) 


where I is the nxn identity matrix. 


(2.4) can be simplified further as follows: 


-1 

H = (I + G^G ) G G 
pc pc 


(I + G„G^)'^ (I + G„G^) - (I + G„G )'^ 
pc p c' pc 


I - (I + G,,G ) 
'■ pc'' 


-1 


(2.S) 


(2.5) shows that the elements of H depend in a very simple way 

on the cofactors of the elements of the matrix I + G G , i.e,. 

pc* * 


(I+G G^) . . 
h, . « 1 - - P 

det(I+G„G ) 
p c 


i — 1 , * , n 


( 2 . 6 ) 


(I + G G ) . . 

h. . « - J P_c^Ji. 

13 det (I+G^G ) 

p c 


i» j = 1, * * * ,n 

i ^ i 


(2.7) 


2ihe argument s will be dropped whenever no confusion exists. 
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where det(I + G G^) denotes the determinant of the nxn matrix 

p c 

I (I '*■ GpG^)^^ denotes the cofactor of the jith clement 

of I + G G . 

pc 


Let the inputs to the system be polynominal inputs with only one 
term, e.g., step, ramp or parabola inputs, which are of primary 
importance. The Laplace transform of each input Tj (t) . j - 1, 
• * • , n, is then 


r, (s) = L (r . (t) ) 
j j 



( 2 . 8 ) 


where r (without argument ) is a constant, and k. is a positive 

integer, e.g., if the jth input is a step then k.=l, a ramp then 

J 

kj=2, etc. 

Then, by (2.7), the steady-state decoupling criterion (2.3) 
becomes : 


1 im 
s-*-o 


n 

z 

j-1 


s- 


s ^ 


-(I + G^G ) . . 

det(I+G G ) 
p c'^ 


= o 


for all i=l. 


,n 


Thus, the following fundamental theorem for steady-state decoupling 
is developed: 
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theorem 2.1: 

Assume that the given plant G (nxm) is stabilizable through the 

p 

configuration of Figure 2.1, then the system is steady^state 
decoupled if and only if: 

o for all i»l , * ’ • ,n 

(2.9) 

Theorem 2.1 can be simplified further for Systems whose inputs 
are not fixed. This is desirable in mo$t practical applications, 
for example, consider an aircraft as our plant G , the thrust, 
flap and elevator inputs must not be fixed in order to perform 
different functions. 


det(I.GpGc) 


Thus, for inputs with arbitrary constants r., 


we have 


THEOREM 2.2: 


Assume that the given plant G^Cuxin) is stabilizable through the 
configuration of Figure 2.1, and that the constants Tj in all the 

inputs Tj(t), j=l'**,n are arbitrary, then the system is steady- 
state decoupled if and only if 


9 



0 


for all and if^j 


( 2 . 10 ) 


PROOF: 


a) Necessity: 

Suppose there exists i^,j' such that (2.10) is not true, then, 
by choosing r^..(t) as the only non-zero input, we have 


n 

1 im z 
s^o j-1 




det(I.GpG^) 




1 im 
s-^o 




sO'j'-iJ 


J 


det(I+G G ) 
P c 


( 2 . 12 ) 


By our hypothesis, (2.12) is non-zero, hence (2.11) is non- 
zero, then hy Theorem 2.1, the system is not steady-state 
decoupled. 


b) Sufficiency: 

Since (2.9) is simply a sum of (2,10) for different values of 
i j, if (2.10) is true, (2.9) is obviously true, hence the 
proof. 


Q.E.D. 


3 

with the jth input, both Theorem 2.1 and Theorem 2.2 can be 
applied. to systems whose inputs are either all of the same type 
(e.g. all inputs are steps) or hybrid (e.g. input 1 is step, 
input 2 is ramp, input 3 is parabolic, etc.) 

Both Theorem 2.1 and 2.2 are in neat mathematical fotms. However, 
they cannot be applied directly, since our objective is to 
determine specifically what to put in the matrix G as the 
compensator functions in order to decouple the steady states of th 
system. Therefore, further result than (2.9) and (2.10) is nec- 
essary. 

Direct approach, which utilizes the expansions of both the 
determinant and the cofactors of a matrix, is used. A general 
result will be given in Chapter 4. Before going into the general 
problem, however, two simple cases are treated first in the 
following chapter. 
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3« SIMPLE CASES 

In this chapter, 2-input, 2-output and 3-input, 3-output plants, 
both compensated by diagonal using the feedback configuration 
Figure 2.1, will be considered. 

Details for the 2x2 case are presented in Section 3.1. Then, in 
Section 3.2, the outline and results for the 3x3 case are given, 

3,1 2x2 CASE 

For a given 2-input, 2-output plant, 

®pll ®pl2 

G = 

®p21 ®p22 

if the diagonal compensator matrix 
^cll 0 

G = 

® 0 ®c22 

is used, the system configuration in Figure 2.1 becomes 
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By Theorem 2.2, for arbitrary 


r^(t) and 

r2(t). 

the system is 

only if 



liin 

1 


s**-o 

s(kj-i) 

det(I+G G ) 

p C-' 

and 

1 

. ^'*Vc>21 

s-^o 

3^2-1) 

<iet(I+G G ) 
P c 


constants in both of the inputs 
steady-state decoupled if and 

(3.4) 

(3.5) 


where k and k are 


defined as in (2.8). 


Let the inputs r^(t) and r 2 (t) be two step functions with arbitrary 
amplitudes, then by (2.8), - 1 and r^.Tj are two arbitrary 

constants. 

Then, by substituting (3.1), (3.2), and (3.3) into (3.4) and (3.S), 
we have 


lim ^p2lScll 

s-*o f^*Spjjg^jj + gp22gc22*®pll®p22Scli8c22 ~ 


®pl2®p21®cll®c22^ 
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and 


lira . 

s-o (l*gpjig 


Note that in both 


^Pl2^c22 

cll^®p22^c22*^pll^p22^cll^c22 

*’^pl2S21^cll^c22^ 

(5.6) and (3.7), the s factor 



from the 


(3.7) 


Final 


Value Theorem was cancelled by the i factor in the input transforms, 
hence no explicit powers of s appears in (3.6) and (3,7). 


Thus, for systems as in Figure 3.1 with arbitrary amplitude step 
inputs, the necessary and sufficient conditions for steady-state 
decoupling are (3.6) and (3,7), 

For a given plant, all the are known, hence the design for 

steady-state decoupling is simply the determination of g , , and 

c 1 1 

^c22* that (3.6) and (3.7) are satisfied. 

For example, consider 

Gp(s) = 


1 1 

s+1 s+2 

1 

s + 3 


By (3.6) and (3.7), if 


1 

s + 4 


15 


g (s) 

cll 


^c22 






(3.8) 


where g" (s) and g' ^^(s) do not contain any pole or zero at the 
cll c22 ' ^ 

origin, or alternatively, lim g' and lim g' ape non-zero 

s->-o <=“ s -^0 *=22 

finite constants, then 


gcijCs) .gc22(s) “ 


c22 


(s) 


Since, 


III (*pllS22'®pl2®p2l)= 1 


n 0 


(3.9) 


the following terra in the denominator of both (3.6) and (3.7), 
^®pll®p22®cU®c22 ■ ®pl2S21®cll«c22^ contains a l/s^ 

factor, will go to infinity faster than both of the numerators in 

(3.6) and (3.7) as s approaches zero. Thus, (3.6) and (3.7) are 

satisfied and the system is steady-state decoupled. It is seen 

that the pure integrators in and g and the constraint 

Cll C 2 2 

(3.9) are important. These constitute the highlights of the 
analysis that follows. 
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Let 


g‘ . = s g" 

Pij PiJ 


g . . = s g' . . 

* ClJ 


(3.10) 


where t . t are integers that will he referred to as the 
P^j cij 

type numbers of the corresponding transfer functions, and 


^ nii’ ® that: lim g' and lim g' .. are non^ 2 ;ero 

^ s->o pij s->o 

constants (i.e. the numerators and denominators of g' . and g" 

Pij cij 

do not contain powers of s as their factors). Whenever g .. =0 
^cij- corresponding K"pij and are defined to 

identically zero, however and become indefinite in this 

case, and we will use the symbol x to identify them for reasons 
that will be clear in Chapter 4. 


The matrices Tp = = O^cij^mxn referre4 

to as the TYPE NUMBER MATRICES of the plant and compensator 
respectively. 


For example, given 

5 

I 

I 

s 

s + 2 


2 

s (s+3) 
0 
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the type number matrix is 


T « 
P 


0 

-1 


By separating the powers of s in each of the transfer functions as 
in (3.10), (3.6) and (3.7) can be expressed as: 


lim 

s-^o 


lim 

s^o 


’ ^^p21*^cll^ 


* ^ p21 • ^ c 1 1 


f pl2*^ c22 


= 0 


= 0 


(3.11) 


(3.12) 


where 




8 pll« cU*® 


-(‘p22"*c22> 


8'p22E'c22 


pll p22 cll c22-'„. , , 

8 pUE p22^ Cll® c22 

_^-(*pl2*'^p21**cll^*c22\. . 

® P12® p 21 ® cll® c22 (3.13) 


Thus, the necessary and sufficient conditions (3,6) and (3.7) 
assume different forms in (3.11) and (3.12). Again, once a plant 
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is given, g ... and hence t » g' . . are known. Therefore, only 

pij PiJ 

and g'cij (5*11) and (3,12) are left adjustable. It was 
shown in a previous example that pure integrators in and 

gc 22 (see (3.8)) are important in steady-state decoupling. In 
terms of the expressions given in (5.10), this is the same as 
saying that the values of 3^^*^ ^c22 factors in the 

attainment of a steady-state decoupling scheme. 


In. order to find out the constraints on t . . and t such that 

cll c22 

(5.11) and (3.12) are satisfied, the following theorem is developed 


THEOREM 3.1 


Let 


a) 


CqCs) 


k 

i = is 


-n 


Pi(s) 


i -m . 

l+rs ^q.(s) 
j=l J 


(3.14D 


be a rational function in Si where P.(s), q.(s) are themselves 

1 1 


rational functions such that lim P.(s) and lim q . (s) are non^zero 

s-»^o s-^o ^ 

finite constants and n^, are integers for all i = li*’‘ik. 
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b) N4Max 


M4,Max {nij , ‘ } 


k 

c) lira I P.(s) ^ 0 (3.15) 

s-^o i = l 

n . =N 
1 

I 

lira I q j (s ) ?^ 0 
$■♦■0 j = l 
m =M 

j 

Z 

lira t (s) -1 

s -^o j * 1 

lDj=0 

k 

where J) Pj^(s) means that the summation is only over those P ■ in 

isrl 1 

n.*N 

1 

(3.14) that have s”^ as their multiplication factor. The other two 
summations in (3.16) and (3.17) are defined similarly. 


Then lim C(s) = 0 if and only if 
s-»-o 


either M>N 


(3.18) 


(3.16) 


(5.17) 


or 


M<0, N<0, 
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PROOF : 


Cq (s) can be expressed as 


Since 


k -n . k ~n . k -n . 

Is ^P.(s) + rs P.(s) + Is ^P.(s) 
1 ^ ^ 1 J- 


c (s) 

o 


i = l 
n . >0 

-j 


i = l 

n . =0 



i^l 

n^<o 


£ *in . 


£ -m . 


£ «m . 


l+Es ^q.(s) + Es ^qj(s) + Es ^q.(s) 

j=i ’ j=i ^ j=i j 


>0 


nij = 0 


nij <0 


k -n 


£ -m 


lim Es ^P,(s)=lim Es ^qj(s) 
s-^o i-l ^ s-»-o j = I ^ 


= 0 by (a) 


n . <0 
1 


m . <0 
J 


lim C (s) 
s-^o ® 


k -n^ k 

lim Es p, (s) + lim j;Pi(s) 

s*^o i = l ^ s^o i = l 

n^>o n^*o 


l + l im is q , (s) 
s->o j = l ^ 

m . > o 
1 


+ 


£ 

lim iq.[s) 
$-►0 


The limit value can be determined for each of the following nine 
poss ib 1 e cases : 
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1 ) 


N>0, M>0 

k M-n. k 

lim IS ^P^(s) + lim is”p.(s) 
s-»>o i = l s-^o i*l ^ 


1 im 
s-^o 


C 

o 


(S)- 


1 im 
s-»-o 


n . >0 
1 


n . =0 
1 


I M-m. 

s^+ lim Is ^q.(s) 
S'^'O j-1 ^ 

m . >0 
J 


I 

+ lim Ss^q . (s) 
. §-*-0 js»l ^ 

m . »o 
1 


k M-n. 

lim IS ^P.(s) 

$■*’0 i-1 
n.>o 

± = 0 if and only if M>N, by (3.V6) 

lim J:q.(s) 5 (3. IS) 

s*^o j=l 

m.=M 

J 


2) N>0, M«0 

k 

lim s’^ IP , (s) 
s->o i=l^ 

lim C (s)= 

s-»-o ^ i 

1 + z<\^ (s) 
j=l^ 
m . =0 
1 


by (3.15) 


3) N>0, M<0 


lim s”^ IP.(s) 
s->-o • i = l^ 

lim C (s)_ ^i''^ 

S-»-0 O - 1 


:s +C0 


by (3.15) 
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4) N=0, M>0 

k 

1 im Z P . (s) 

s-^o i = l ^ 

n . =0 
1 

lim C (s)= = 0, 

s->o ^ , . -M- / \ 

1 im s Zq . (s) 

m . =M 
3 


by (3.16) 


5) N=0, M=0 


k 

lim ZP.(s) 
s-^o 1*1 
n . *0 

lim C (s)= ^ 

s*^o ® ^ 

1+lim Zq.(s) 

s-^o j = 1^ 

m . =o 
3 


/ 0, 


by (3.15) 


6) N*0, M<0 

k 

lim ZP.(s) 
s-^o i = l 

lim C fsl _■ "i°° 

® 1 


^ 0, 


by (3.15) 


7) N<0, M>0 


lim C 
s->-o 


o 


(s) 


0 


M 

lim Z q. (s) 

S-^O js=l ^ 

m.=M 

3 


0 , 


by (3,16) 
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8) N<0, M*0 


1 im 
s-*-o 


C 

o 



1 4-1 im Eq , (s) 
s-^o j = 1 j 



0 , 


by (3,17) 


9) N<0, M<0 


lim C (s)_ 0 

s->-o ^ 1 


Thus, lim C (s) = 0 if and only if one of (1), (4), (7), (8), 
s-^o ^ 

(9) is true. Since the conditions in (1) to (8) are equivalent 

to lim C (s) = 0 if and only if M>N, and (9) gives M<0, N-;0, the 
s-^o 

theorem is proved. 

Q.E.D. 


Note that (3.18) contains only the powers of s, neither p. nor 

appears in this expression. Also, note that C^(s) is of exactly 
the same form as the rational functions in (3.11) and (3.12), 
therefore, the theorem can be applied directly. 
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Compare (3.11) and (3.12) with (3.14), and by the definition of 
M and N in (b) of Theorem 3.1, we have: 

I 

Si2*%21*"c11*S22’ 

N_=Max {t .. +t } * t +t^ii 
12 p21 cll p21 cll ^ 

Njj^Max (tpi2*tc22^ ' *pl2*^22 


Where the notation Max {*•’} denotes the maximum value among all 
the elements in the brackets and the subscripts on N are used in 
accordance with (3.2) and (3.3) to distinguish them from each other. 


Since known for a given plant, the only unknowns in 

(3.19) are and which can be chosen to satisfy (3.18). 

Once t ^22 chosen, M, Nj 2 » ^21 known, and (3,15), 

(3.16), (3,17) can then be written down explicitly. 


In general, these expressions contain both g" .. and g^ . Since 

pi 1 ci j 

the pole and zero locations and the gains of each g^ are free 

c 1 j 

parameters, they can hopefully be adjusted to satisfy (3,15), 
(3.16), and (3.17). 
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These free parameters should also be designed for stability and 

transient response of the system, therefore, they cannot be 

adjusted with complete freedom. However, as was mentioned before, 

(3.18) does not depend on g' therefore, the design of stability 

c 1 j 

will not destroy the steady-state decoupling as long as (3.15), 

(3.16) and (3.17) are not violated. Hence, once t .. and t are 

c 1 1 c 2 2 

determined, stabilitycanbeconsidered. 

After all are designed, however, (3.15), (3,16) and (3.17) 

must be checked. If satisfied, the design is completed, if not, 
slight adjustments of the free parameters, under the allowance of 
stability, can be made in order to satisfy these constraints and 

i 

hence guarantees that the steady states are decoupled. 


It might happen that in some cases, no adjustment in . . is 

cij 

possible to satisfy these constraints, e.g.. 


1 im 


S'+'O pll^ p22“S pl2^ p21^ 8 cll8 c22 ^ ^ possible if 


sio ^^"pll®"p22'^pl2^p21^ * ^ happens to be true for the 
given plant. 


In cases like this, another choice of t^j| and t^22 is necessary 
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Following this procedure, we can, at present, assume that (3.15), 
(3.16) and (3.17) are satisfied. Then, by Theorem 3.1, if the 
constraints (3.15), (3.16) and (3.17) for the rational function in 
(3.11) are satisfied, then (3.11) is true if and only if 

either ^^^12 

(3.20) 

or M<0, ^<0 

12 

Similarly, if the constraints (3.15), (3.16) and (3.17) for the 
function in (3.12) are satisfied, then (3.12) is true if and only 

if 

either M>N 2 j 

(3.21) 

or M<0, N2 j<0 

For steady-state decoupling, both (3.11) and (3,12) must be true, 
therefore, combining (3.20) and (3 . 21 ), each of the following four 
setsofcriteriacanbeused: 


M>N 

M>N 


12 

21 


(3.2 2) 
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M>N 


12 


M<0, N2j<0 


(3.23) 


M>N 

M<0 


21 


N 


12 


<0 


(3.24) 


M<0 




(3.25) 


It should be noticed that (3.23) and (3.24) are redundant since 
they are contained in (3.25), 


The best choice among these four sets will depend on the type 
number matrix of the given plant. 


Consider the following example: Given the 2x2 plant^ 


G (s) 
P 


-s +1 -s +2 

(sTTj^ (s+i)2 


- 3s •♦•1 “ s + 1 

3(s+1)2 Cs+1)2 


The plant is taken from an example in (5) . 


(3.26) 
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Find a steady-state decoupling scheme using diagonal G and 

c 

the configuration of Figure 2,1 assuming that the inputs are 
arbitrary steps. 


By inspection, tp.j = 0 for all i, j == 1 , 2, then by (3,19), 
M=:Max {t^^j,t^22*^cll'^^c22'^cll*^c22^ 

^12= '^cll 
'^21° ^c22 


If (3.22) is used, fcH = t^22”^ simplest solution (note that 

the solution is not unique). For this particular choice, 5^11* 

8 ^c22~T ^^c 22' “the introduction of pure integrators 

in the loops will cause the steady states to be decoupled. 

(3.23), (3.24) and (3.25) can also be used, however, in this case, 

the solutions for both and ^^22 turn out to be negative, 

which corresponds to the introduction of differentiators in G and 

c * 

is physically undesirable. 
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Since only one term appears in the numerators of (3.11) and (3.12), 
(3.15) is satisfied automatically by definition of 8‘pij 
g’cii C.see (3.10)). 


By inspection of (3.13) and by noting that both of the la$t two 
. -M 

terms contain s for the above choice of t -- and t we have 

c 1 1 c 22 

for (3.16) 


lim 

s*^o (g pllS ’ p22^ * ell® ' c22 “ ^’pl2^*p21^'cll^*c22^ ^ ^ 


Since 

lim , . . 

$-►0 pll^ p22 " ^'pl2^'p21^ ■ ^ ^ 

We have 

lim , , / 

s-^o ®cll^c22^^ 

which is again satisfied automatically. 

Similarly, by inspection of (3.13), (3.17) is also satisfied auto- 
matically, since by the above choice of t^jj and t^22^ none of the 
terms in (3.13) has 0 as the power of the associated s factor. 
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Therefore, we are guaranteed to have a steady-state decoupled 
system by introducing one pure integrator in each of and 

®c22' 

Actually, in this case we don’t need (3,15) and (3,17), since 
M=2>0, Nj 2®N22=1>0, and by (1) in the proof of Theorem 3,1, 
only (3,16) is sufficient. 

It should also be noticed from the proof of Theorem 3,1 that the 
constraint (3.15) was used only to make (3,18) also a necessary 
condition. If (3,15) is not true, the sufficient part of the 
theorem is still guaranteed by (3,16) and (3,17), Therefore, 
it is usually only necessary to check (3,16) and (3.17) in 
practical design. 

For inputs other than steps. Theorem 3.1 must be generalized as 
follows; 

THEOREM 3.2 

Let 

C^(s) = 1 - Cp(s) (3.28J 

s 

where t * 0, 1, 2, Cq(s) is as defined in (3.14) and (b) , 

(c) are the same as in Theorem 3,1. 
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\ 

( 

(3.29) 


the result follows immediately by Theorem 3.1, 

Q.E.D. 

Now let input 2 be a ramp, while input 1 is still a step (i.e., 
kj=l, k2=2). Then by (3.4), (3.6) and (3.11) are still the same. 
However, for input 2, since k2=2 in (3.5), there will be an 
additional s factor in the denominators of (3.7) and (3.12), 

Since the only difference is this additional s, (3.19) remains 
unchanged, and the application of Theorem 3.2 gives M>N 2 j+l. 
Therefore, the conditions corresponding to (3.22) becomes: 


Then 


1 im 

S**‘0 ^ t ^ ^ ^ 


= 0 if and only if 


either M<0, N+t<0 


or 


M>N + t 


PROOF: By writing: 


CJs) = 


Is 
- i*l 


k -(n^ + t) 


Pi(s) 


t -ra . 

1 + Zs Jq. (s) 
j = l J 
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M>N 

12 

For the plant (3.26), by (3.27) we have t , =2, t =1 as the 

cll c22 

simplest solution. 

Thus, we need one more integrator in in order to decouple the 

steady states, if input 2 is a ramp instead of a step. 

3.2 3x3 CASE 

For a given 3-input, 3-output plant, 

®P11 ®pl2 ®P13 
®p21 ®p22 ®p23 

^p31 Sp32 ^p33 

if the diagonal compensator matrix 

8cll ° ° 

0 8c22 0 

° ° «c33 

is used, the system configuration is as shown in Figure 3.2. 








FIGURE 3.2 


Expressions for det(I+G G ) and (I+G G ) , , can be obtained either 

pc p c ij 

by direct expansion as was done in Section 3.1 (see (3.1), (3.2) 
and (3.3)), or by using the formulae (4.5) and (4.6), then by (2.10) 
of Theorem 2.2 and assuming step inputs (i.e. k =k =k,*l),a set 
of 6 limit expressions similar to (3.11) and (3'.12) can be 
obtained. Compare these with (3.14), we have 

M ^%ll*^cll’*p22*^c22’^p33'^^c33’*^pll**p22*^cll*^c22> 

*^pl2'^*^p21*'^cll'^^c22>^pll*’^p33*^cll*^c33>^pl3'^’^p31 
*‘cll‘"*^c33>^p22*^p33*tc22**c33>’^p23*^p32'"^c22*’^c33> 
^pll*^p22*^p33**^cll*'^c2 2*'^c33’*pl2*^p23**p31*^cll 
^*c22**c33’^pl3**p21*^p32*’^cll''*^c2 2'''^c33''^pl3''^p22 
^p3l'^^cll*'^c22*'^c33’^pl2**p2l'^*p33**cll**^c22'*^*c33* 
^pll*^p23*^p32*^cl 1*^022*^^033^ 


(3.30) 
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"12 = ”^=' ^%21^%ll’*p21*%33^^n*"c33' 
^p23*’^p31*^cll*^c33^ 


13 p31 cll p31 p22 cll c22* 


^p32'^^p21*^cll'^^c22^ 


N2j=Max ^tpi2*^c22'^pl2'^^p33'^^c22*^c33* 
^pl3'^^p32'^^c22*^c33^ 


N-_=Max {t ,^+t *^,t ,^+t ,,+t 
23 p32 c22 p32 pll cll c22 

t ^ +t ^ +t +t' } 

p31 pl2 cll c22 


(3.30) 


N_ =Max {t +t +t ^ +t +t 

51 pl3 c33' pi3 p22 C22 c33* 

^p23'*’^pl2‘*’^c22'^^c33^ 


N, =Max {t +t ,t +t +t +t 

p23 c33 p23 pH cll c33 

+t - 1 +t } 
p21 pi 3 cll c33 


Note that ^ for some i,j in a given plant, any term in 

(3.50) that contains the corresponding t has to be drooned 

paj 

In Chapter 4, an analytical scheme will be designed to take care 
of this. 


s 
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Similar to (3.22) through (3.25), we have 64 possible 

sets of criteria here to choose from. However, similar to the 
previous case, only the two corresponding to (3.22) and (3.25) 
are not redundant , these ar e : 


M>N 


M>N 


12 


13 


M>N 


M>N 


M>N 


21 


23 


31 


M>N 


32 


and 


M<0 


N 


N 


12 


13 


N 


21 


N 


23 


N 


31 


N 


32 


<0 

<0 

<0 

<0 

<0 

<0 


(3.31) 


(3.32) 


Consider the following example : 
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EXAMPLE: 

Given a 3x3 plant^ with 

g (s) « 0.081(s-0.20SHs + 0.967 + jl.379) 
(s+0.967-jl.379)/D*D 

g IjCs) = -6.12(s + 0.837Hs + 0.947 + jl.l44) 
(s+0.947-jl.l44)/D-D_ 
g 13 ( 8 ) » -202(3+1.885) (s-13.037)/D-D^ 

g 2 i(si “ -0.00163(3+2.881) (s+0.032+j0. 313) 

(s+0.032-j0.313)/D-D 

gp 2 2 (s) • -0. 153(3+0. 8 24) (s- 0.047 + jO. 205) 

(s-0.047-j0. 205)/D-Dp 
gp 23 Cs) = -9.07(3+26.339) (s+0.03+j0. 361) 
(s+0.03-j0.361) /D-Dg 
gp 3 l(s) = -0.00209(3-1.049) (s + 0 . 268) / 0- D 
gp 32 (s) = 0.0995(3-0. 12) (s+3. 485) /0-Dp 
gp 33 (s) = -235. 5(3+0. 361+jO. 076) 
(s+0.361-j0.076)/D-Dg 


(3.33) 


NASA STOL C- 8 A aircraft, with thrust, flap angle and elevator 
angle as the inputs and velocity, angle of attack, pitch angle 
as the outputs. 


57 


where 

D(s) .= (s + 0, 018 + j 0.336) (s + 0 , 018 - j 0 . 336) 

(s+1.103+j 1.277) (s+1 .103-jl. 277) 

D (s) = (s+0.99+j0.479) (s+0 . 99- j 0 . 47$) 

TH 

0^(s) =* (s + 3.3 + jl0.49) (s + 3,3-jl0.49) 

E 

Dp(s) = s + 1 

Find a steady-state decoupling scheme using diagonal G^, and the 
configuration of Figure 3.2, assuming that the inputs are 
arbitrary steps. 


By inspection of (3.33), t .. =0 for all i i = 1,2,3. Thus (3.30) 

pij ' 

assumes the following simple form: 


M=Max {tcll’*c22’*c33’’^cll**c22’’^clX*'^c33’ 
^c22**c33’*c11*^c22'^*g33^ 

Njj=Max ^tcii.tcxi+tc22^ 

Nji'Max {tc22'*c22'^*c33^ 

N23=Max I^c22‘^c22*^ctl^ 

Njj=Max <tc33.tc33+tc22> 

Nj2=Max Uc33,t^33+tcii} 


(3.34) 
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By (3.31) and (3.34), it is clear that M must be 1 1 *^c22^^c^53 * 

and tcll> ^c22’ *^c33 (otherwise t^,j j + + 

cannot be a maximum). Hence the simplest solution is * 

tc 22 “ ^c33 " means that the introduction of one pure 

integrator in each of the compensators 8c22* ®c33 cause 

the steady states of the system to be decoupled , 

Since M>0, N^^>0, only (3.16) has to be checked. It can readily be 
found that this is satisfied, therefore we are guaranteed to have 
a steady-state decoupled system. 

(3.32) can also be used, however, as in the 2x2 example of 

Section 3.1, the result requires pure differentiators in 

gc22 Sc33* bence also not desirable for this particular plant. 
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4. GENERAL nxm CASE 

Results of Chapter 3 show that, for the two special cases 
considered, the s t eady-state , decoupl ing criterion can be written 
as a set of n(n-l) inequalities, where n is the number of outputs 
of the plant (or number of inputs or outputs of the system). 

These results will be generalized in this chapter to systems 
consisting of m-input, n-output plant G^(nxm), n-input, m*»output 
compensator G (mxn), and unity feedbacks are employed as shown in 
Figure 2 . 1 . Exactly the same approach as in Chapter 3 is presumed 
and it will be seen that both Theorem 3.1 and 3.2 are applicable. 

As was shown in Chapter 3, the first step is to obtain expressions 
for det(I+G G.) and (I+G G^).. as in (3.1), (3.2) and (3,3). For 

jp ^ jP ^ J 

the general case, this can be accomplished by using the following 
formulae which are proved in Appendix A. 

For any nxn square matrix G 

n 

det (I+G) =1+1 I 

Jl = l If ki < * * * <k„ 

1 I 

n-2 

(I+G).. = - I r 

^ Jl = 0 lskj^< • • • 

' • * ,k^j^i, j 


ki , • • ‘ ,k 
in \ k, , • • • ,k 


(4.1) 


G i,k. , • • • ,k 



(4.2) 
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/ ii,* ’ * ,io \ .... 

where G| I denotes the minor formed from rows 


and columns matrix G with 1 < £.< n. 


Let G - GpG(. where Gp is the nxm plant matrix and G^ is the mxn 
compensator matrix. Then by Binet-Cauchy formula (6), we have 


k j j ‘ ‘ > k 

ki , ‘ ‘ • ,k 


ki , * ‘ • ,k 




kj. • • • .kjj 


i.kj, • • • ,k 

j » k i , j k 


Z G 

l^p < • • • <p ^ P 
0 


1 * kj^ , • • • , 

Pq>Pi» * * * > P 


pQ»Pl» ' ' * »P 
j » k 1 > i k . 


)t<m- 1 


i>m-l 
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Combining (4.1) and (4.3), (4,2) and (4.4), respectively, we 

have 


min (n ,m) 

det(I+G G )=1+ I ^ I 

P ^ Jl = l 1 1 < • * * <k <n li Gi < * * • <a„ 


^ ®1 * j ^ \ ^1 » * » 


(4.5) 


min (n- 2 ,m- 1 ) 


I I 

likj<' • *<kj^«n l^p^<* • »<p^cm 


*>k. 

G I ^ IG 





PQ»Pj, * ’ ’ »PjJ 


(4.6) 


Then, by expanding the associated minors , (4,5) can be written as 


min (n ,m) 

det(I+G G )=1+ Z 


P c 


£*1 l<k. < * • • <k . sn l<oi < • • • <a*sm 
1 t ^ i 


, CTfl /K 

G 


i ^ r i „ /*'A 


o* o'. , • • * ,0 


‘ p 




p a' 


<’1 ’ 


> o 


I 

i a 


'i 



.G 
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where 



a” 

I 


{ 1 if is an even permutation of 

-1 if is an odd permutation of o*' , * ’ • , o" 

* St ^ St 


and E , E i^epresent summations over all possible permutations of 
0 ^ 0 *' 

and respectively . 


The 


identity 6 




used, which can be proved by first rearranging into 

^l* ‘*^i* using simple reasoning. 


It can easily be shown by letting n=m=2 in (4.7) that (3.1) can be 
obtained through this expression. Similarly, general expressions 
for the cofactors + can be obtained by (4.6). Thus, by 

p j 1 * 

using (4.1), (4.2), (4.3), and (4.4), the problem of expressing 

det(I + GpG^) and terms of the transfer functions 

g_. . and g^. . explicitly is solved. 

Pi; cij 
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Then, following the approach of Chapter 3, limit expressions 
similar to (3.6) and (3.7) can be obtained. Ip order to evaluate 
the values of these limits, it was found convenient to express 
each transfer function as in (3.10). By doing so, (4.7) can be 
written as : 


min(n,m) 

det(I + GG)-l+ S X 

P ^ £ = 1 l<ki<’*’<k in 1 soi < * * ‘ <a „ sm 

1 9 1 £ 


I 

a ’ 




Note that there are 


min (n ,m) 

^ n^£ *m^£ • ^ 

£=sl 


£ 


min (n ,m) 

^ n^£ *m^£ 
£sl 


(4.9) 


terms in (4.8). Similarly, (4.6) can be manipulated into the 
following form: 
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(1+ G G ) . 
P 


min (n-2 ,m- 1) 

E . I I 

£=0 l<ki<*--<k <n 1 spq<‘ • 

kj. • • • j 


• I z 

p 'p" 





C*. 10) 


Again , note that there are 
min (n - 2 , m- 1 ) 

min(n-2,m-l) 

= E *n-2^£ + l 


terms in (4.10). 


(4.11) 


By (4.8), (4.10) and Theorem 2.2 (assuming arbitrary T^), n(n-l) 
limit expressions similar to (3.11) and (3. 12) can be obtained. 

These limits are, according to Theorem 2.2, necessary and sufficient 
conditions for steady-state decoupling. In order to satisfy these 
conditions. Theorem 3.1 and 3.2 were developed to find the 

cij- 


constraints on t 
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By comparing (4.8) and (4.10) to the denominator and the numerator 

of (3.14), it can be seen that they are of exactly the same form, 

only that n^,m^,p^,q^ assume more complicated forms here. 

Therefore, similar to what was done in Chapter 3, constraints on 

t .. for steady-state decoupling can be obtained by applying 
c 1 j 

Theorem 3.1 (or Theorem 3.2, if the inputs contain ramps or parabolas 
besides steps) to each of these n(n-l) limits. 


To be more precise, let’s go through these step by step as 
foil ows ; 


1. Assume step inputs with arbitrary ampl itudes , (i . e . , kj*!, r^. 
arbitrary for all j-l,’*‘n), and consider the configuration 
Figure 2.1. By Theorem 2.2, the system is steady-state 
decoupled if and only if 


lim 

det(I*GpG^) 


for all i,j=sl,**‘,n and i^j . 


(4.12) 


Note that 


- 1) i 1 since k. 


for step 


inputs . 
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2. By comparing the denominator and the numerator of (3.14) to 
those of (4.12) which are given respectively in (4,8) and 

(4.10) , we have J m’s, J q's, L n's and L p's (see (4.9) and 

(4.11) ) as follows: 



where each possible combination of k’s o's and p's under the res 
trictions in (4.8) and (4.10) contributes to one of the above. 


3. Then, for each of the n(n-l) limits (4.12), Theorem 3.1 can be 
applied and a set of inequalities consisting of 


Max (m , I j '=1, • • • ,j} 


and Max (n.,| i'=l,''-,L) 


(4.17) 


(4.18) 
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can be obtained as in Chapter 3. Again, since all the T’s in 

(4.13) and (4.14) are known for any given plant, the only 

unknowns are the T's which can be chosen to satisfy the 

c 

inequalities and hence achieve a steady-state decoupling scheme 


4. Whenever any transfer function in or is identically zero, 

those terms in the summations of (4.8) and (4.10) that contain 

such a factor will also be identically zero, hence the number 

of non-trivial terms in (4.8) and (4.10) will be less than J and 

L, respectively. The number of m and n^, will also be 

reduced. Thus, those m. , and n , in (4.17) and (4.18) associated 

J ^ 

with the identical zero term should be dropped, since they don't 
even appear in (3.14). In order to express this analytically, 
the identification symbol "x” introduced in Chapter 3 (see the 
discussion following (3.10)) will be used. Also, the following 
definition of annihilation sum is needed: 


The annihilation sum is defined to be a summation, which will sum 
up to be an empty set whenever there exists at least one 
identification s ymb ol x in the summands , otherwise it is the 
same as algebraic sum. The symbol (+) will be used for such 
kind of summation, e.g. 

l( + )2( + )3 1 + 2 + 3 = 6 

l(+)x(+)3 = 0, an empty set. 
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By using these concepts, (4.13), (4.14) become 



Now, starting from the type number matrices T and T (see 

p c 

Chapter 3), we know immediately from (4.19) and (4.20) which 
m^ , and n are to be included and which should be discarded. 

5. Define W, N. . to be the IxJ and IxL vectors with their elements 

J ** 

corresponding to the J and L annihilation sums given in (4.19) 
and (4.20) ♦ Note that some of their elements can be an empty 
set, whenever is contained in those particular elements. 

In this way, whether a term should be dropped or not is expressed 
analytically. Then , by Theorem 2.2 and Theorem 3 . 2 , the follow- 
ing general theorem for steady-state decoupling can be givep: 

Theorem 4: Let the given nxm plant be compensated by an mxn 

as in Figure 2.1. T^ (nxm), T^ (mxn) are the type number 
matrices of the plant and compensator, respectively. M, N.. 
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are defined to be the maximum among the elements of M and . , 

J i 

respectively. Where M and ^ are the IxJ, and IxL vectors 
defined above, then under the following constraints: 


(i) 


min(n-2,m-l) 

lim I I I I 

s-^o i = 0 lsk,<’»*<k„sn l<p^<*‘*<p sm p» 

<N..> ^ ^ 0 ^ 

kj, ,k£?^i,j 



( 4 . 21 ) 


min (n ,m) 



min (n ,m) 
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Where <N. > under the summation sign of (4,21) denotes that the 
Ji -N,. 

sum is only over all those terms with s as their 

multiplication factor in (4.10). The <M> and <0> in (4,22) and 
(4.23) are defined similarly. 

The system is steady-state decoupled if and only if either 

for all i = j 

Theorem 4 looks formidable, however, foi^ systems with less than 
4 inputs and 4 outputs, long-hand calculation is still feasible, 
especially when assumes some simple forms like being diagonal, 
as were shown in Chapter 3, which is often of practical 
importance . 

Besides, due to its analytical nature, Theorem 4 can be programmed 
into computers, thus making the design easier. 


M>N. . + k.-l 


or M<0 N. < 1-k. 

J i J 





PART II 


STABILITY AND DESIGN 



51 


5. INTRODUCTION TO PART II 

Part I gives the scheme for decoupling the steady states of a 
system, however, it should be noticed that: 

1. The result does not guarantee stability, 

2. The whole discussion is meaningful only when the 
closed -loop system is stable. 

Therefore, stability must be considered after the steady-state 
decoupling scheme is achieved. 

The problem of stability and design of multivariable systems has 
been widely investigated (e.g. (5), (7), (8), (9), (10), (11)), 

and a survey of the existing methods was given by Anderson (12). 
In general, efforts have been made to utilize the beauty of the 
existing single-loop techniques such as Nyquist-Bode-Nichol ‘ s 
methodsandrootlocusdesign. 

In this part of the thesis, a new connection between single-loop 
and multivariable systems is seen by properly factorizing the 
closed-loop characteristic equation. This makes the design of 
multivariable systems possible by using any suitable classical 
single-loop method. An extended root locus method is developed 
and diagonal will be considered primarily. 
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6, DESIGN OF CLOSED-LOOP SYSTEMS WITH 
2x2 PLANT AND DIAGONAL G 

c 

The characteristic equation for single-loop systems with cascade 
coBipen sa t ion and unity feedback is 

1 + g„(s)*g^(s) ^ 0 ' (6.1) 

r 

where g(s) is the given plant and g(s) is the cascade compensator 

if c 

function to be designed. 

Two major techniques for the design of g (s) are Bode's method and 

c 

the root locus method (see (14)). However, only one unknown 
function can be handled in each of these methods. Therefore, they 
cannot be applied directly to multivariable cases, since in general, 
there exists n*m unknown compensator functions to be designed. 

A simple factorization which shows the connection between single- 
loop and multivariable cases will be given in this chapter. It c^n 
then be seen that the above-mentioned single-loop methods can still 
be applied for multivariable systems. 

The design philosophy will be illustrated through an example in 
Section 6.4. Before that, however, several important steps must be 
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established. These are given in Section 6.1 to 6.3 as follows: 
6.1 CHARACTERISTIC EQUATION 

It is proved in (7) and (8) that the stability of a multivariable 
system as shown in Figure 2.1 is determined by the zeros of 


Nj (s) 

A 

and N(s) 

, which are defined as 

follows: 


Nj(s) 

A det(I+G (s) -G^(s)) 



(s) 

(6.2) 


N(s) 

Di(s) 

(6.3) 

Where 

the rational function N^(s)/D^ 

(s) is in irreducible form. 


i.e., no common factor between Nj(s) and D^(s) is left uncancelled. 
And A^,(s) Ap(s) represent the characteristic polynomials of the 
rational transfer function matrices Cs ) and Gp(s), respectively. 

The characteristic polynomial of a proper rational transfer function 
matrix G(s) is defined to be the least common denominator of all the 
minors (in irreducible rational form) of G(s) (see e.g, (13)). 

Therefore, by (4.5), if all the minors of G and G are non. zero 

p c 



54 


and no common pole exists between G and G*. the denominator of 

P ^ 

det(I + GpG^) is A^.(s) Ap(s) before any pole-zero cancellation is 

performed. Since Nj(s)/Dj^(s) is in irreducible form, it is clear 

by (6.3) that N(s) simply consists of all those factors that were 

cancelled in getting the irreducible Nj(s)/D^(s) form. Therefore, 

in this case, if all the common factors in (6.2) are left 

uncancelled, the zeros of (6,2) alone determine the stability of 

the system. Unfortunately, the same conclusion is not true in 

general if zero minor(s) of G and (or) G exists. Furthermore. 

P c 

if ^'all^' the common factors are left uncancelled, erroneous results 
can still be ‘obtained as was shqwn by Chen in (7). 

However, it is found (see Appendix B) that if 

1. cancellations are selected systematically by using ( 4 . 5 ), 

2. poles of are carefully selected, 

then (6.2) alone determines the stability of the system. 

Since (1) above can always be done and (2) can be taken care of 
fairly easily in the process of design, the mathematical 
possibilities in which zeros of (6.3) must be considered can be 
bypassed . 
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Thus, det(I+GpG^) = 0 (6,4) 

will be referred to a:: the characteristic equat i on for multi- 
variable systems as shown in Figure 2.1. 

Details are given in Appendix B , 

6.2 CONNECTION BETWEEN SINGLE-LOOP AND MULTIVARIABLE CASES 

For 2x2 plant G^, and 2x2 diagonal G^, 

detCI + GpG^) = '^p5®cn8c22 

which can be factored as 


det(I.GpG^) = (l*gpllgcll) 

r,. gp22*C<^®«p)gcll 

L l^gpllgcll J 

(6.6) 

^^*®pll®cll^ 


(6.7) 


f 1 


where G , «p 22 ^ Cd^tG^) g, , ^ 

- ®p22 L ^P22 

(6.8) 


^■^^pll^cll 


By (6.7), the roots of the characteristic equation are simply the 
zeros of the rational function 1 2c 1 1 ^ ^c22^ * 
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It is also clear by (6.6) that for non-trivial cases (G 0) , 

eq 

all the zeros of the first factor will be cancelled exactly by 
some of the poles of the second factor. Therefore, the roots of 

■ «c 22 = 0 (6-9) 

will determine the stability of the system. 

The similarity of the form of (6.9) to that of (6.1) suggests 
immediately the possibility of applying the single- loop methods 
mentioned above to multivariable cases. But, unlike the 
gp in (6.1), which is known for a given plant, of (6.9) is not 

a known function . Hence, neither Bode plot nor root locus for 
^eq drawn at this stage. By inspection of (6.8), the only 

unknown function contained in is We can, of course, 

choose an arbitrary function for then the only unknown 

function left to be designed is g^22* design is reduced to 

that of the single-loop case. For example, let be an 

arbitrary constant, say 3, then for any given plant, can be 

obtained through (6,8), and the design of g^22 carried 

through by using (6.9). 

However, in doing so, G may turn out to be very unstable, which 
will make the design of extremely difficult. Therefore, 

instead of choosing it arbitrarily, a guide in designing g^^^ is 
preferable. 
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By inspection of (6.8), it can be seen that the roots of 


1 + 


det G. 


g 


g . - = 0 

^cl 1 


p22 


( 6 . 10 ) 


and 


1 + g 


pll^cll 


ss 0 


( 6 . 11 ) 


constitute part of the zeros and poles of G , respectively. 

e cj 


Again , both • (6 . 1 0) and (6,11) are of the standard form (6.1). 
Furthermore, detC^/g^^^ and are now known functions. 


Therefore, any single-loop method be used in designing 

to place the roots of (6.10) and (6.11) at desirable locations. 
Since these roots will be part of the poles and zeros of G 

eq 

what is meant by placing them at desirable locations is that 


g 


cll 


should be designed such that these roots, together with the 
other known poles and zeros (see (6.8)) form a reasonably good 
pole-zero pattern for G^^ (i.e., G^^ is not badly unstable). 

Once g^jj is designed, all the poles and zeros of G^^ are known, 
and the problem of designing is, by (6.9), reduced to that 

of the single-loop case, and can be done by either Bode^s or 
root locus method. 
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In summary, what has been accomplished so far is that the effect 

of on the pole-zero pattern of and hence on the design 

of g can be seen through (6.10) and (6.11). Therefore, (6.10) 
C Z 4 

and (6.11) serve as a guide in designing g in order to make 

V X A 

easy the design of g . 

C Z 6 


Because of the standard forms involved, both Bode*s and root locus 

1 

design techniques can be applied. For better insight of the 
problem, the root locus method will be considered primarily. 


6.3 . IDENTIFICATION OF POLES. AND ZEROS 


Since the design will be concerned with (6,9), the poles and zeros 
of ^gq(s) must be well identified, and the problem of pole-zero 
cancellation must be considered carefully. This can be done 
generally by considering the sum and ratio of the following two 
rational functions G^^(s) and 0 ^( 5 ) ; 


Gj(s)4 


Nj(s) 

oTTsT 


( 6 . 12 ) 




N,(s) 

5717 


(6.13) 
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Where N^(s), D^(s) are the numerator and the denominator polynomials 
of G^(s) (i=l,2),also note that N^(s), Dj(s) here are different 
from those in (6.2), 

For simplicity, the argument s will be omitted in the following 
discussion. 


Let Dj 2 greatest common factor between and D 2 , and 

^11* ^22 remaining factors as shown below: 

Di = Djj Di2 

C6.:4) 

= ^22 °12 

Then, the sum of G ^ and G 2 is 



= ^1 ^22 ^2 ^11 
Dij Dj 2 D22 


(6.15) 


Note that ^12 ^22 the least common multiplier of Dj 


and D 2 * 
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Now, consider the sum 1 + (6.12), (6,13) and (6,14) 


1 



1 , "2 . 


= 1 + 


- 

^22 ^12 


N, 


(6.16) 


If the common factor 0^2 between Dj and is cancelled, and no 
other cancellation is performed, (6.16) becomes 


1 + ^ - ^22^1 * ^11^2 

Gi D22 Nj (6.17) 

The numerator of (6.17) is exactly that of (6.15). Therefore, 

the zeros of G^^ + G , before any possible cancellation by the poles, 

^2 

would be the same as those of 1 + ,if (and only Di.,) is 

Gj I*-' 

cancelled. No other cancellation in G 2 /GJ is allowed, even if it 
can be done. Otherwise, some zero of (6.15) would not appear in 
(6.17). 

Therefore, if only is cancelled in forming G 2 /GJ, the root 

locus for G 2 /GJ would give all the zeros of G^+G^ if no 
cancellation is done between the numerator and the denominator of 
(6.15). 
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Let = { * * ’ } and Zq = {’*•} denote the set of poles and zeros of 

* 

the rational function G, where the multiplicity of each pole or 
zero is. counted, e.g,, 

let G(s) = 

s(s + 1)2(s + 2) 

H(s) = (s*3)^(s*5)^ 
s2(s+1)3(s+2} 

then 

P^. = {0,-1, -1,-2} 

Z = {-3, -3, -3, -5, -5} 

( 6 . 18 ) 

Pjj = {0,0, -1,-1, -1,-2} 

= {-3, -3,-5, -5} 

Also» let s^ n =2 denote the intersection of the two sets s^ and 
$ 2 , defined as in set theory only that multiplicity is taken into 
account here,e,g.,in(6.18) 

n Ph " (0, -1,-1, -2} 

Zq n Zh ° {-3, -3, -5,-5) 
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Furthermore, s^+s^ is defined to be the set consisting of all 
the elements of s^ and those of $ 2 , (counting multiplicity), e.g., 
in (6.18) 


P 


+ 

G 


H 


{ 0 ,- 1 , - 1 ,- 2 , 0 , 0 , - 1 ,- 1 , - 1 ,- 2 } 


Similarly, Sj^-S 2 is defined to be the set formed by taking away 
all the elements (counting multiplicity) of s^ from s^, and is 
defined only when S 2 is a subset of s^* e.g., in (6,18) 

Ph - Pg = (o.-n 

Zj. - Z„ = {-3} 

By using these notations and by inspection of the denominator of 
(6.15), we have 


•■gi.gz ' (’’GiH Pep 


(6.19) 


Similarly, by comparing the numerators of (6.15) and (6.17), we 
have 


ZG1+G2 = ^1 



+ 


( 6 . 20 ) 
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The poles and 

zeros for the ratio 

G 2 /GJ 

can also be 

obtained 

in a 

similar way by 

inspection 

of the second 

term in (6. 

16) 


PG2/G1 

- Gj G 2 

(Pgi n 

PG3) 



(6.21) 

^G2/Gj 

= ^G2*PGj - 

(’’gj n 




(6.22) 


Note that in (6.19), (6.20), (6.21) and (6.22), D^2 is cancelled, 

and no other cancellation is performed. 

The application of (6.19), (6.20), (6.21) and (6.22) will be 

illustrated through a design example in the next section. At 
present, however, (6.20) will be used to justify one statement 
pointed out in Section 6.2, i.e., the roots of (6.9) will determine 
the stability of the system. 

Let = l+gpiigcii 

then (6.5) becomes 

det(I+GpG^) = Cj+G2g^22 


(6.23) 
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By (6.20), 


‘^1*°2®c22 


Z G 2 
1 + ^ ®c22 


if the common factors between the denominators of and ^2^c22 
(and only these common factors) arecancelled. 


Since G^/G^ 


G by (6.8), 


2, G2 

*^Gj 8c22 



^c22 


Hence, the zeros of (6.23) are exactly the same as those of 

• gj,22» ^ ® roots of (6.9) are the same as those of the 

characteristic equation (6.4). Therefore, they do determine the 
stability of the system (under the restrictions given in 
Appendix B) . No pole would be lost on account of the factorization 
and the using of (6.19), (6.20), (6.21) and (6.22). 

6,4 DESIGN EXAMPLE 

Consider the2x2plant 


G 

P 


s-^3 

s(s+T) 

3 

s + 2 



(6.24) 
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For diagonal the characteristic equation (6.4) is, by (6.5) 

%^®cllgc22 ' ° (6.2S] 

By the factorization (6.7) and the discussions in Section 6,2 and 
Section 6.3, the design philosophy follows: 

1. Design according to (6,10) and (6,11), to achieve 

a reasonably good pole-zero pattern for 

2. Design g^22 ^<^cording to (6.9) to meet system specifica- 
tions. 


where (6.9), (6.10) and (6.11) are repeated below: 


1+G 


eq *^c22 


0 


i detG„ 

1+ E 

S22 ^<^11 


0 


^■^^pll^cll ° 


(6.9) 

( 6 . 10 ) 
( 6 . 11 ) 


By (6.24), 

a s S ^3 

s(s+l) 


(6.26) 
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(6.27) 


d e t G = 
P 


14s^ + 10s-^12 
s ^ (s + 1 ) ( s + 2^ 


(6.28) 


detG 
P 

®P22 


14s^^-10s-»-12 1 ^ 
_ s^(s+l) (s*2) _ 


(- |) 


_ 7s^+Ss*6 
s (s + 1) (s + 2) 


(6.29) 


Note that the common factor s between the denominators 
of det Gp and g^22 cancelled. Also note that the use 
of diagonal is allowed in this example, since 
multiplicity of each plant pole is not reduced in the 
minor det G^ (see Appendix B) . 

(6.30) 

(6.31) 

then, by (6.8), 



‘"1 ° ^'^®pll®cll 
°2 = 


(6.32) 
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The poles and zeros of are by (6.21) and (6.22), 


•"eq ' ^Gj) 


= ^G^+^Gj-cPCj 0^0^) 

The poles and zeros of and G 2 are by (6.15)and (6,20), 


(6.33) 


‘’°l = ’’gpllgcll = C6-34D 


l + g_, ,g 


pll^cll 


(6.35) 


0 

■ ’’gp22*‘’C<l®tGp)g^lj - [ ‘’gp22 ^ ’’f‘^®*‘^p)«cll] 

= {0}+{0,0,-l,-2}.Pg^jj.^{0} n (tO,0,-l,-2}+Pg^jpl 

= {0,0, 0.-1. -2} + Pg^ll-{0} 

= {0.0,-l,-2} + Pgj.i2 (6.36) 


®p22 


.»cll 


(6.37) 


By (6.34) and (6.36) , 


pg. n P 


G, = (O.-D + Pg 


cll 


- P, 


(6.38) 
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2}+P„.i-K0,-l}+P 1 

L ScllJ 

C6.39) 
(6.40) 

Thus, has two known poles at' s-0 and s = -2, the other poles are 

the roots of 1+gpx i^cl l“^ * which is (6.11). The zeros of G are 

simply the roots of 1 ^ ^‘Scll * which is (6,10). 

^p22 

Then, the design procedure fol lows : 

1. Prepare the root loci for 

- detG^ 

and 1 + ki^ £ =0 

gp22 

where kj is a real parameter. The result is shown in Figure 6,1 
and 6.2. The loci for positive and negative values of kj are 
represented by solid and dashed curves, respectively. 


Then by (6.33), 


’’g ' g 

eq ^pll^cll 


Z, +{0,-2} 

^gpl l^c 1 1 


'p22 


^ g_,, ®cii 
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2. The roots of Figure 6.1 will be poles of G , therefore, for 

eq 

negative values of k-, there will be one pole of G in the 

1 ^ eq 

right half plane (on branch©). Similarly, by Figure 6.2, 
a zero of will be in the right half planeCon branch ©) if 

kj<0. This is certainly undesirable. Thus, negative kj will 
not be considered. 

3. Mark down the known poles s=0, and s=-2 of G on Figure 6.3 

eq 

and superimpose the loci of Figure 6.1 and 6.2 corresponding 
to positive kj on top of it. Note that roots on branch © 
and @ correspond to poles and zeros of respectively. 

4. Increase the value of k^^ from 0 to «, and observe the change 
of pole-zero pattern. It can be seen that 

(i) 0<kj<<5 is not desirable, since the poles will be 

clustered together near the origin. 

■ (ii) kj>>5 is also not desirable, since the pole on branch® 
will be pushed into the negative real a^^is, hence, 
dominant roots of (6,9) will probably be determined by 
the two known poles s=0, s*-2 and the zeros on branch ©, 
which are too close to the origin. 
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5. When the roots on branch © which are the zeros of 

will be close to the two zeros on Figure 6.2 and the pole on 
' branch © is as shown. Another zero, by Figure 6.2 will be on 
the negative real axis at about s = -35. This pole-zero 
configuration looks to be the best, since it is possible to 
confine the roots of ~ corresponding to the two 

poles s =0 and s =-2 and the two conjugate zeros, to be on the 
negative real axis . And at the same time , therootof 
l'*'^ 2 ^eq " ® branch starting from the pole at kj = 5 will 

be someplace to the upper left of the pole, which is a good 
location for dominant root. Therefore, try k^^sS, 


6 . Once kj=5 is determined, all the poles and zeros of are 

known . The root locus gain (see Appendix C) of G^^ , denoted 
by k is found as follows: 

^ eq 



-2(s-Zi) (S-Z2) (3-23) 

(s-pi) (S-P 2 ) (S-P 3 ) (S-P 4 ) (6.41) 

i 
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where , and are the zeros and poles. By inspection of 
(6,41) and the definition of root locus gain, we have 

K == -2 
eq 

7. The root locus for l+k^G^ =0 is then drawn as shown in 

2 eq 

Figure 6.4. The choice of is now strictly that of a 

single-loop problem. It is easily found that the two small 

real roots meet each other at about k^k^ = 0,36, i.e. k.,» 

eq c 

0.36 ^ , 

k * For this value of k«, the smallest real root 

eq I 

on branch (J) wil 1 be 'at approximately the breakaway point and 
will hopefully be the best among all the possible locations 
on this particular branch. The root on branch @ corresponding 
to J^2^eq ' 0.36 is also shown in Figure 6.4, It is seen 

that this is a pretty good dominant root. 





73 


The schematic for the designed system is as shown in Figure 6.5. 



FIGURE 6.5 

The simulation result is shown in Figure 6.6, in which r^^^lO and 
T2=5 were used as reference inputs and the time responses for the 
two outputs Xj and are as shown. 

It is clear that stability has been achieved. The transient 
response of X 2 is very good, however, that of y^ is kind of slow. 

If this is not allowed by the specification, a redesign is necessary. 
However, as in the single-loop case, every trial, despite of its 
failure, provides some guide for the next trial. 
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In this example, it can be seen from Figure 6.4 that the roots on 
branch 1 are, compared with that on branch 2, very close to the 
origin,. This is probably the cause of the slow response exhibited 
in y^ , What we can do is put pole(s) and zero(s) in (instead 

of just using a pure gain) to push the roots of Figure 6,2 farther 
away from j u axis. Then the two conjugate zeros in Figure 6.4 
will also be farther away from the jw axis and the resultant root 


locus for 


G 

eq 


will move toward the left, 


transient response. 


thus improving the 


Incidentally, the steady states, of the outputs are decoupled. This 
is because gpjj and gp22 ^ proved very 

easily by (3.19) and (3.22). 
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7. DESIGN OF CLOSED-LOOP SYSTEMS WITH 
3x3 PLANT AND DIAGONAL 

7.1 GENERAL 

It will be shown in this chapter that the same design philosophy 
given in Chapter 6 can be carried over to 3- input , 3 -out put 
plant. Again, the central idea is in the factorization of the 
characteristic equation det C I + GpG^.) * 0 . 

For 3x3 G and diagonal 3x3 G ,,the system configuration is as 
shown in Figure 3.2 and the characteristic equation is by (4.5), 

det(I + GpG^) = ^'*'2pllgcll'^®p22®c22'^®p33®c33‘^^^p^332cll^c22 

^^p^ 22^cll^c33*^^p^ ll^c2 2^c33*^*'‘^^^^p^®cllSc22Sc33 
* ® (7.1) 

where (Gp)jj, (Gp) 22 » denote the cofactors of gpU* gp 22 

and gp 33 > respectively. 



77 


(7»1) can be factored as 


fl4g. g 1 ri^ gp22^C<^p)33gcii ^ gp33MGp)22g cll 

1 + p --P *c22 


L ^'^Spiigcii 


l+gpllgcll 


(G )ii+(detG )gcii ^ „ 

“Vn ^ Sc22gc33 = 0 

i*gpllgcll J 


By defining. 


‘^l^^*®pll8cll 
°2^p22*^'’p^33gcll 
®3^p33'^ *-°p^ 22gcll 

(7.1) and (7.2) become 


^I'^^2®c22^^3gc33'^^4®c228c33 * ^ 


fi.^i 

L Gj 


1 = 0 

®c22 g,®c33 Gj®c22 *c33 


Sc33 

(7.2) 

(7.33 

(7.43 
(7. S3 
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Similar to the 2 x 2 case considered in Chapter 6, the roots of 
( 7 . 4 ) will be exactly the same as those of: 

G G G 

^•^T^8c22*r^8c33'^-^gc22ec33 ' ° 

Gi Gj G2 

if the common factors between the denominators of G2 and Gj, 

Gj and Gj, and (and only these common factors) are cancelled 
in G2/G2, G3/GJ and G^/Gj, respectively. 

Thus, if the poles and zeros of G2/G2, ^3/6^ and G4/G2 are obtained 
through (6.21) and (6.22), which were designed to meet the above 
cancellation restrictions, the roots of ( 7 . 6 ) are the same as those 
of the characteristic equation (7.1). 

By defining 

F I - ^P 22 *^°ph 3 gell 

^ °1 ^*8piigcll 

r. ^3 ®p 33 "^ 22®cll 

F, A — - 

“ ^*8piigc„ 

^ ^ ^1 " ^"«pll8cll 


( 7 . 7 ) 


( 7 . 8 ) 


( 7 . 9 ) 
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(7.6) becomes 


**’’i®c22*’'22c33*''38c22Sc33 ° ° 


(7.10) 


(7.10) is of the same form as (6.5), therefore, the same 
factorization can be done on (7.10) to give 


(1-Flge22) • 



^2*^38c22 

^*'’l®c22 


®c33 


= 0 


(7.11) 


Again, by (6.20), if the common factor between the denominators of 
^2^^32 c 22 and l'*‘Pigc22 (^^^ only this common factor) is cancelled in 


G A ^ 2^^ 5^c22 

“ 1*Pi8c22 ^ 


(7.12) 


the roots of (7.10) will be the same as those of 


1+G *g,,=0 

eq c33 


(7,13) 


Since (7.10) is simply (7.6), and the roots of (7.6) are the same 
as those of the characteristic equation (7.1), if (6.21) and 
(6.22) are used in determining the poles and zeros of Fj, and 
Fj, the following conclusion can be made: 



80 


The roots of (7.13) are the same as those of the characteristic 
equation (7.1), if (6.19), (6.20), (6.21) and (6.22) are used in 

determining the poles and zeros at each stage. 


Thus , stability design can be considered through (7.13), The 
form of (7.13) is exactly that of a single-loop characteristic 
equation (6.1). Therefore, similar to the design of Z^22 
2x2 case in Chapter 6, any single-loop design method can be applied 
in designing once all the poles, zeros and the root locus 

gain of G__ are known. 


The ingredients of the poles and zeros of G can be found by 

. eq ^ 

substituting (7,7), (7.8) and (7,9) into (7.12) to express G in 

eq 

terms of the elements of G and G explicitly. However, due to 

P . 

the fact that some cancellations must be done while some others 
are not allowed, this approach may sometimes lead to an erroneous 
result. Therefore, the analytical schemes (6.19), (6.20), (6.21) 
and (6.22) designed to take care of the pole-zero cancellations, 
are recommended. 


By applying (6.21) and (6.22) on (7.12), we have 

’’‘^eq ° ' ^’’’=2*^38022 ’’l*Flgc22> 

z z 

^eq ” ^2'*’^3^c22 ^"^^l^c22 ^^^2‘^^3®c22 ^ ^^'^^l2c22^ (7.15) 
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(7.14) tells us that the roots of 



(7.16) 


(7.17) 


constitute part of the zeros of G . 

eq 

These justify what can be seen by inspection of (7.12). Actually, 
these were done by inspection in Section 6.2 (see (6.8), (6.9), 

(6.10) and (6.11)), before the development of (6.19), (6.20), 

(6.21) and (6.22). 
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Thus, root loci for (7.16) and (7.17) identify part of the poles 
and zeros of and similar to the 2x2 case in Chapter 6 , 

these root loci can be used in the design of (note: this 

corresponds to Chapter 6 , compare (6.10), (6.11) with 

(7. 16) and (7.17)). 

However, unlike the 2x2 case, these two root loci cannot be drawn 
directly, since Fj and F 3 /F 2 depend on which is also to be 

designed (hence not known yet!). The dependences of Fj, F 2 and 
^3 ^cll given in (7.7), (7.8) and (7.9). 

Again, by repeated application of (6.19), (6.20), (6.21) and 

(6.22) on (7.7), (7,8) and (7.9), or simply by inspection, it can 

beseenthattheroot loci for 

^ * ®pll 8 cll =0 C7.183 

and 1 + - ^3^ g . 0 

gp 22 (7.19) 

give part of the poles and the zeros of F , respectively. 
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And the root loci for 


1 * -^gcll = 0 

®p33 


(7.20) 


and 1 + P e ^ = 0 

cl 1 


(G ) 


11 


(7,21) 


give part of the poles and the zeros of F /F«, respectively, 

3 * 

Now, let^s look back and see what we've got: 

1. We concluded that roots of (7.13) are the same as those of the 
characteristic equation (7.1). Therefore, the design of g 


c33 


can be done through (7.13) if G is known. 

eq 


Some of the poles and zeros of G are adjustable through g 

eq ^ *c 22 ' 

and the relationships are given in (7.16) and (7.17). Thus, 

the effect of g on the pole-zero pattern of G and, hence, 

eq 

on the design of 8 ^. 33 , can be seen through (7.16) and (7.17). 
Therefore, (7.16) and (7.17) serve as a guide in designing 

®c 22 design of not formidable. This is exactly 

what was obtained in Section 6.2 for the 2x2 case. 
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3. Since both (7.16) and (7.17) are of the standard single-loop 

form ( 6 . 1 ), 8^.22 designed if Fi* and P 3 /F 2 known. 

4. Again, some poles and zeros of Fj and ^ 3/^2 adjustable 

through and the relationships are given in (7.18), 

(7,19), (7.20) and (7,21). Thus, the effect of on the 

pole-zero pattern of Fj and F 3 /F 2 , and hence on the design of 
^c22* through (7,18), (7.19), (7.20) and (7.21). 

Therefore, these four equations can be used as a guide in 
designing 

Thus, it is clear that (7 . 16) , (7 . 17) , (7.18), (7.19), (7.20) 
and (7,21) are important in stability design. By the standard 
forms they assume and by their similarity to those in Chapter 6 , 
it can be concluded that root loci for these six equa t ions can 
help us design g^jj and g ^»22 ^ stable enough G such that 

gc 33 designed according to (7.13), 

Since gpj^» ^p22* ®p33' ^^p^ll' ^^p^22* ^^p^33 detGp are known 

for a given plant, the four uncompensated root loci (i.e., g i^sk, 

a free parameter) for (7.18), (7.19), (7.20) and (7.21) can be 
constructed right after the plant is given. However, the other two 
root loci for (7.16) and (7.17) cannot be drawn until after g^j. 
is designed, since both Fj and Fj/Fj depend on g^jj. As mentioned 



85 


before, and ^ 2 some poles and zeros other than those 

given by (7.18), (7.19), (7.20) and (7.21). Therefore, the 
identification of these poles and zeros is necessary, both for 
constructing root loci for (7.16) and (7.17) and for guiding the 
design of This constitutes the topic of the fol lowing 

section. 

7.2 IDENTIFICATION OF POLES AND ZEROS 


As mentioned in the previous section, direct algebraic manipulation 
may lead to an erroneous result., so let's apply the analytical 
schemes (6.19), (6.20), (6.21) and (6.22) to identify all the poles 

and zeros weareinterestedin. 

The expressions for all the poles and zeros of have already 

been given in (7.14) and (7.15). The poles and zeros in the sets 
77 

1 + Figc 22 ^2‘*'^3^c22 roots of (7.16) and (7.17), 

respectively, and can be taken care of by the corresponding 
root loci. The other poles and zeros left to be identified are 
respectively the elements of the following two sets: 


(7.22) 
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In order to express these explicitly in terms of the poles and 
zeros in the plant and the compensator repeated applications 

of (6.19), (6.20), (6.21), and (6.22) on 

?!, F^, F 3 /F 2 necessary.. The results can be written down by 

inspection as follows: 


^ 2 *^ 3 ^ c 22 ° '''’2 * * *’^022 ’ ^'”''2 G ’’>'36022^ 


'’^*'’l®c22 ' ’’^1 " '’®c22 


''1 ®p22**^'"p^33®cll 

(■ • ■ 

{’’l + gpllgcll ’’gp22'^('^p^33®cll^ 

° ^8p22"(Gp)338cll * '’^^SpllScll 
' ^'’l*8pll8cll ^’’Sp22*^S^33Scll> 
**2 ^'*’®pll®cll ®p33* ^*^p^ 22®cll 

^ ^*®pll®cll ^ ®p33* ^*^p^ 22®cl 1 ^ 


(7.24) 
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^2 ®p33* 22 ®c 11 ^*®pll®cll 

^ ^*®pH®cll ^ ®p33*^'^p^22®cll^ 

(7.24) 

■ <'’'*8pii8cn 

I 

^’’3 ° * '’^*8pll8cll 

■ ^•’I'gpllgcll ^'’( = p)ll*C<l®«p)gcll> 

where 

^p22 ^ 33^c 1 1 ^p22 ^^p^33 ®cll 

^ ®p22 ^ ** ^*^P^ 33®cll ^ 

®p33*^*'p^22Scll ®p33 ^^®p^22 ®cll 

(7.25) 

^ ®p33 ^ '”f°p^222cU^ 

’’(Vll*fdetGp)g^jj = P(Gp)j/P(detGp)-Pg^jj 

■(^Gp)j^n P(detGp)g^ll) 
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P +P 
Spll ^cll 


Also 

^Fj/Fj = ^F3"Pf2'c'’f2 ^ '’"s^ 


(7.25) 


(7.26) 


By proper substitutions of (7,25) into (7.24), then (7.26), (7.23) 
and (7.22), the poles and zeros of interest can be identified. 


This will be illustrated through a numerical example in the 
follow in gsecti on. 


At present, however, some simplifications on the above general 
expressions can be made. It is observed that if. 


P 


P 


P 


g C ‘’(G ) 
P22 P 33 


g ^ ^ (G ) . 
®p33 P 22 




(7.27) 
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where A Cl B denotes that set A is a subset of set 
plicities of the elements in each set are counted , 


p n P = P 

*p22 ^‘^p^33®cll ®p22 

n '■(detGp)g^„ = P(Gp)ii 

®p33 ^ ^‘’P^22 ®c 11 Sp33 

then (7.25) becomes 

'"S22*^°p^33®cll * '"^°p^33*'’®cU 

P P ^ P 

®p33*^‘’p^22®cU ^'"p^22 8^11 

‘’CGp)„*(detGp)g^ll ° 

P- = P +P 

^*«pii«cii «pii ^Cii 


, again multi- 
Then, we have 


(7.28) 


(7.29) 
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and 

’’l + gpllg<-n ’’gp22^^'^p^33®cll ' '’gcll*^’’®pn ^’’('^p^33^ 

'’.1 + gpllgcll ^ *”gp33* ('^p^ 228 cll ° **8cll ^ Spll ^ f*^p^22 

'’l + gpjlgcii n '’(Gp)ii+(detGp)gcu = '’gcii*^'’gpll ’’CdetGpJ 

By (7.29) and (7.30), (7.24), (7.26) become 

'’Fl = ^l*Spllgcll‘"’(S^33'^'’gpll n ’’(Gp)33> 

^’’l = ^8p22*^°p^338cii*'*®pll ^*'®pll ^ /^‘^p^33 

’’F2 = ^l + gpllgcll*^('"p)22 ^''Spll ^ '”^^p^22^ 

^F2'^8p33* V 228 cll*'’ 8 pll’^’’®Pll 

’’F3 = ^l*gpllgcll*^f‘’"*''p)‘^’’«pn ^ 

= ^ '’(detGp)> 


(7.30) 


(7.31) 
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and 

-tPg^^^n P(detGp) ^'’(detGpr(’’gpiilT’CdetGp))}n 

(7.52) 

■^%il^’’f"p^22^’ ^”(0etGp)-('’gp,2nP(detGp))>n 

^’’CGp522'(^pll^'’(S’22)^ 

Once a plant is given, ^ known set. 

Therefore, by of (7.51), all the poles of Fj are well 

identified. They consist of all the roots of (7.18) which are 
adjustable through some other fixed poles given by the 

above set which is known. 

Similarly, by inspection of (7.31) and (7,32), all the zeros of 
Fj, all the poles and zeros of F 3 /F 2 are well identified. Some 
of them are fixed , the others , which are the rootsof(7.19), 
(7.20) and (7.21) respectively, are adjustable through g^jj* 
Therefore, can be designed according to the four root loci 

for (7,18), (7.19), (7,20) and (7.21), to realize a presumably 
good pole-zero pattern for Fj and ^ 3 /^ 2 * such that the design of 
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^c22 (7.16) and (7.17) will not be formidable. 

As anyone who is familiar with root locus design knows, there is a 
certain amount of tr ial-and-error involved. This is more so in the 
multivariable case because of the successive dependence of the root 
loci described so far. However, a little experience can always 
lead to good judgements that would reduce the amount of the trial- 
and-error. For example, in the design of described above, if 

the root loci for (7,18), (7.19), (7.20) and (7,21) extend well 
into the right -hal f -plane, more sections for g^. j j is , in general , 
recommended. Otherwise, it is very probable that the resulting 
^1 ^3^^2 both) contain poles and zeros well in the right- 

half-plane, thus making 8^.22 <^ifficult to design. This is, of 
course, a trade-off between 8^.22* more sections of 

compensation are used in ^3^^2 be made more 

stable, hence, less sections are required in 8^.22* Conversely , if 
gcii is chosen to be too simple, more sections should be needed in 
gc22* Judicious choice can be made by investigating the four 
root loci (7.18), (7,19), (7,20) and (7.21), and can be 

designed accordingly. 


The poles and zeros of G can also be identified in a similar 

eq 

way. This will be clearer after the consideration of the numerical 
example in the fol lowing sec t ion . 
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7.3 DESIGN EXAMPLE 

As was found in Section 7.1 and 7.2, g^.j2 designed according 

to the four root loci 


1 + C g =5 0 

^pll^cll 


(7.18) 
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where ^ given in (7.7), (7.8), (7.9), and 

(7.12), respectively. 

The general design procedure then follows: 


1 . 


Identify all the poles and zeros of 


F 3 /F 2 and 


2. Prepare root loci for (7.18), (7.19), (7.20) and (7.21) with 

gcii * free parameter. 

3. By varying the value of kj from to ®, observe the 
accompanying changes in root locations . 


4. Choose the value of that corresponds to the best pole-zero 
pattern for V and F 3 /F 2 sense that 8^22 designed 

most easily to give good pole-zero pattern for G . 


5, If no value of kj gives satisfactory and F 3 /F 2 , use pole- 
zero pair as necessary in g^^^ to pull the loci toward the left 
and determine the gain value for best pole-zero pattern for 
Fj^ and ^ 3 /^ 2 * 


6 , Construct the root loci for (7.16) and (7,17), usine b « =k-. 

&c22 2 


a free parameter. 
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7. 


Adjust 1^2 in 3 to find best pole-zero pattern for 
in poles and zeros as necessary, as in S. 



Put 


8. Construct root locus for (7,13). 


9* Design to meet specifications. 


Consider the 3-input, 3-output plant (3.33) repeated below: 

gpU (s) =0.081 (s -0.205) (s+0.967+jl .379) 

(s^0.967-j 1.379) /D*D^ 

t 

gpjjCs) =-6.12(5+0.837) (5+6.947+) 1.144) 

(s+0.947-jl.l44)/D-Dp 
8pl3(s)=-202(5+1.88S) (5-13.037)/D'Dg 
■gp2j(s) =-0.00163(5+2. 881) (5+0.032+) 0.313) 

(s+O.O32-)O.313)/D-03,jj 

*p22^®^” -*153(5 + 0.824) (s-0.047+)0. 205) (7 

(s-0.047-)0.-20S)/D-Dp 
g 23 (s) =-9-07 (5+26. 339) (5+0.03+) 0.361) 

(s+0.03-)0.361)/D-Dg 
gp 3 l(s)=- 0 . 00209(5-1. 049) (5+0. 268) /D-D 
gp32Cs) =0.0995 (5-0. 12) (5+3.485) /D*D 
g (s)=-235.5(s+0.361+)0.076) 

if O i 

(s+0.361-)0.076) /D-D_ 

c 
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where 


D(s)=(s+0.018+j0.336)(s+0.018-j0.336) 
(S+1.103+J1.277) (s+ 1.1 03- j 1.277) 
D (s)=(s+0. 99+ j 0.479) (s+O . 99- j 0 . 479) 
Dg(s)=(s+3.3+j 10.49) (s+3. 3- jlO.49) 
Dp(s)=s+1 

By manipulation. 


,36.85(5-0.096) 


%^ll°gp22gp33-gp23gp32°-- W^6 


,-19.08(5+0.229) 


^°p’22 = gpll®p33-gpl3gp31=^T^S 


E “TH 


ffi 1 „ ,-0. 0224(5^+1. 656S+0. 694) 

^V33 gpllgp22 gpl2gp21 ^ 


detG = 5.232 

P D-Dg.Dp.D^^ 


By inspection of (7.33) and (7.35), it is clear that ( 
is true. Therefore, by (7.31), 


’’'=1 = 'l*8pll8cir'°F 


''l ' gp22*^‘^p^338cll 


(7.34) 


(7.35) 


. 27) 


(7.36) 
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^2 ° ^®p33*^''p^22®cll 


1+g g *^Db*2d_ 
pll CH E F 


^^3 = 


where ^D = { (-0 . 018+j 0 . 336) , (- 1 . 103±j 1 . 277) } 
^Dth = {-0-99±j0-479} 

Zp = {-3'.3±jl0.49} 

\ = {-!) 


are the sets of the zeros of D (s) , D_ (s), Dp(s), D-,(s) 
1 TH r 

(7,34), respectively, 


and note that 



2 ,^ + 2 ,. +2 
F 


D D'ri ^TH 


etc . , 


and by (7.32) or by (7.26), (7.37) and (7.38) 


^F3/F2 = 'CGp)3i,CdetGp)g^jj 


(7.37) 


(7.38) 


in 


(7.39) 
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Then, by (7,23) 

' + g It8 1 1 75 

2 3^c22 ^pll^cll E F *c22 

Finally, by (7.14), (7.15) and (7.40), we have 


= ^1+Fig *^D. 
eq 1®C22 E 


= ^F,+F,g 
eq 2 3^c22 


( 7 . 40 ) 


( 7 . 41 ) 


Theoutlineof thedesign thenfollows: 

1. Prepare the four root loci for (7.18), (7 . 19) , (7 . 20) and 
(7.21), using These are shown in Figure 7.1 to 

Figure 7.4 (note that dashed curves are the loci corresponding 
to negative kj) . 

By inspection of these loci> the following can be observed: 

(i) kj<0 is undesirable, since for negative kj, there will 

be one branch in each of these four plots that extends 
along positive real axis to + « and this will tend to 
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produce a pole-zero pair on positive real axis for 

G . which is certainly undesirable, 
eq 

(ii) When k^SSO, there will be a pole-zero pair of close 

to the point (0.5, 1.5). The reason is that the roots 
corresponding to ki = 50 on branch (J) of Figure 7.1 and 
Figure 7.2 are close to each other. By (7.36), these 
two roots will be one pole and one zero of Fj 
respectively. Since they are close to each other, the 
root for 1+Figc22^^ corresponding to this pole- zero 
pair is very difficult to push far away from this 
region. Thus, by (7,41), a pole of will be 

around (0.5, 1.5). Similarly, by Figure 7.3 and 
Figure 7.4, there will be a zero of close to the 

same point. Thus, a pole-zero pair of exists near 

the point (0.5, 1.5) in the right-half-plane. This 
will most probably cause the corresponding root of 
(7.13) close to the same point, hence, undesirable. 

(iii) For k^>50, the situation is obviously worse. Therefore, 
the range of that remains to be investigated is 
0<kj<50. 

(iv) It can be seen that when k^ is too close to 0, the root 
on branch 0 of each plot will all be close to the 
origin , hence also not desirable . 



lOQ 


2. When = the pole-zero locations seem to be the best. 


Thus , 


try k^=10. 


3. Using construct the root loci for (7.16) and (7,17) 

with g a free parameter. The results are shown in 

Figure 7.5 and Figure 7.6, respectively. By inspection of 
these two plots, the following can be observed. 

(i) ^^2^^ undesirable for the same reason as that in 
1 (i) above . 

(ii) 0 >k 2>>-10 looks better than the other range, since the 
root on branch in both Figure 7 . 5 and Figure 7,6 
will be farther away from the jw axis, hence, more 
stable G can be expected (note that roots in 
Figure 7.5 and Figure 7.6 give poles and zeros of 
respectively, see (7.41)). 

(iii) Although the root on branch © of Figure 7.5 will be 
close to the origin for 0 >k 2 >>- 10 , it can still be 
tolerated because the same thing does not happen 

in Figure 7.6. Thus, the root of (7.13) corresponding 
to this pole-zero pair of G can still be adjusted to 
be not too close to the origin. 

€ 

Thus, try g^22 * ^^2 " 
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4. The roots corresponding to 8^,22 ” read off from 

Figure 7.5 and Figure 7.6. By (7.14), these constitute part of 
the poles and zeros of Together with the other known 

poles (roots of D^(s)= 0, by (7.41)), the root locus for (7.13) 
can be constructed as shown in Figure 7,7, with 8^22 * ^3* ® 
free parameter. By inspection, = ^3 * ^ value. 

Thus, we have determined all three compensator functions with only 
gain adjustments. The resulting system schematic is as shown in 
Figure 7.8. For the three step inputs r^^ = 126.7 ft, /sec., 
r = -0.25 radian and r, = -0.5 radian, the simulation results for 
the three outputs yj (t) , y2 (t) and y3 (t) are given in Figure 7.9, 
Figure 7.10, and Figure 7.11, respectively. It is clear that the 
resulting system is stable. 







FIGURE 


103 


STOL 

C-8A 

AIRCRAFT 






FIGURE 7.9 


r~ 

1 .^ 


TIME RESPONSE OF OUTPUT IN FIGURE 7.8 WITH 

rj=126.7*u(t) 

T2= -0 . 25 • u (t ) 
r -0 . 5 • u (t ) 


u(t) : UNIT STEP 


FIGURE 7.10 


TIME RESPONSE OF OUTPUT FIGURE 7.8 WITH 

rj=126.7-u(t) 

T2 = -0 . 25 * u ( t ) 
rj=-0.5*uCt) 
uft) : UNIT STEP 
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FIGURE 7.11 

TIME RESPONSE OF OUTPUT y IN FIGURE 7.8 WITH 

o 

r^*126. 7 'u(t) 

. 25 ' u (t) 
r . 5 'u Ct ) 
u(t); UNIT STEP 




PART III 


APPLICATION TO STOL AIRCRAFT 
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8. COMPENSATION DESIGN FOR STOL C-8A AIRCRAFT 
WITH STEADY-STATE DECOUPLING 

8.1 GENERAL 

The simulation results. Figure 7.9, Figure 7.10 and Figure 7.11, 

for the system in Figure 7.8, justify that the root locus 
technique developed in Chapter 6 and Chapter 7, can be used for 
designing both the stability and the transient response of a 
multivariable system. Stability and transient response are 
certainly the most important factors to consider when designing a 
system, however, some other factors are also important. Among 
them (e.g., steady-state accuracy, integrity, sensitivity, etc.), 
steady-state accuracy is usually the most important. In single- 
loop theory, the restriction on steady-state accuracy usually makes 
it impossible to adjust some parameters with complete freedom. In 
root locus terminology, root relocation zones (14) exist, which 
limits some of our abilities to relocate those roots of interest. 

For the multivariable case, due to the existence of mutual coupling, 
the problem of steady-state accuracy is more complex. This can be 
seen by comparing the steady-state values in Figure 7.9, Figure 7.10 
and Figure 7.11 to the input commands . For example , input is 
a step of magnitude 126.7 ft. /sec., while the velocity output y^ is 
only 93.8 ft. /sec. when steady-state is reached. 
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One way to reduce or eliminate the steady-state errors in multi- 
variable systems is decoupling the steady states of the system. 

The concept of steady-state decoupling was developed in Part I 
of this thesis, and the schemes obtained can be used directly to 
determine what types of functions should be used in the compensator. 
Then, results of Part II are applied to stabilize the system. 


8.2 DESIGN PROCEDURE 


The plant under consideration is the longitudinal mode of the 
NASA STOL (Short Take Off and Landing) C-8A Buffalo Aircraft which 
IS a 3-input, 3-output plant as shown below in block diagram form. 


«1 

"2 

“3 


Cthrust, %) 


.flap angle, 
'“radian ^ 


.elevator angle, v 
'“radian ^ 


STOL 

C-8A 

AIRCRAFT 




^2 



.velocity, 
^ft,/sec. ^ 

^angle of attack 
^radian 


^pitch angle, ^ 
radian ^ 


) 


figure 8.1 


Thus, the plant can be represented by a 3x3 matrix G . The transfer 
functions g are given in (3.33) and (7,33). 
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By using the feedback configuration as that in Figure 3.2, a 
steady-state decoupling scheme was obtained in Section 3.2 to be 

^cll"^c22"^c33“^'^‘®*» 

gc22Cs) ' i (8.1) 

Sc33^®^ “ 7 ®c33^®^ ^ 

where, as defined in (3* 10), numerators and denominators of 

^cll (s), ^c22^^^ ^c33^®^ contain power(s) of s as 

their factors . 

(8.1) tells us that the introduction of one pure integrator in each 
of the three compensator transfer functions will cause the steady 
states of the system to be decoupled. Thus, one pole (at the 
origin) is required in each of the three unknown functions 8j.ij(s), 
^c22^^^ ®c33^^^' other poles and zeros and the gain values 

are unknown and have to be determined for stability and transient 
response. The design of these unknowns can be done in exactly the 
same manner as that in Chapter 7, except that the existence of 
the extra pole (s-0) in each of g^.n(s), 8^22^5), and g^^jjCs) has 
to be taken into account. Also note that the pole-zero expressions 
(7.41) (7.36) (7.39) are still valid, since the same plant is 
considered. Then the design procedure follows: 



no 


1) Prepare the 4 root loci for’ (7.18). (7.19), (7.20) and (7.21) 


q k 

with L , where k. is a free parameter. 

ells 


These are shown in Figure 8,2 to Figure 8,5 (note that these 
loci can be drawn by adding the additional pole at the origin 
in Figure 7,1 to Figure 7,4). 

By inspection, it is observed that 

(i) kj<0 is not desirable , for the same reason as 1 (i) of 
the design in Section 7:3. 

(ii) For kj^>0, branch on all four plots are quite unstable, 
which is the effect of the introduced pole at the origin. 
By the same argument as that in 1 (ii) of the design in 
Section 7,3, it can be concluded that no range of k^ is 
desirable . 


Therefore, only one pole (at the origin) and no zero in is 

most probably not enough to give satisfactory system perform- 
ance, some other pole(s) and (or) zero(s) are recommended. 


^The argument s will be dropped whenever no confusion exists. 
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2) By spirule check, or by noting that any introduced zero (that 
goes with the decoupling pole at the origin to form one section 
of filter) on the negative real axis cannot overcome the effect 
of the pole at the origin, it can be seen that if only one 
section of filter (pole at origin, zero on negative real axis) 
is used, the situations will always be worse than those in 
Figure 7.1 to Figure 7.4. Since that design was only marginally 
successful and the situations now are worse, several sections of 
filter are recommended . 

3) By spirule, or by inspection (if experienced enough) it can be 
seen that two sections of filters are enough to pull branch ® 
in each of these four plots into the left-half-plane. Then 

no poles and zeros of and F 3 /F 2 are in the right-half-plane, 
(Note that by (7,36), roots in Figure 8.3 are the zeros of F^. 
Those in Figure 8.2, together with the known pole at s=-l, which 
is the root of Dp(s)= 0 , give all the poles of Fj. Similarly, 
by (7.39), roots in Figure 8.4, Figure 8.5 and the known pole at 
s=-l give all the poles and zeros of F 2 /F 2 ), and the design of 
^c22 ^c 33 continued in the same manner as that in 

Section 7.3. However, by noting that each new root locus that 
has to be constructed out of these four has also an additional 
pole at the origin (comes from g ^^2 ^c33^ ( 8 . 1 )). This 

will tend to destabilize the results and make the design of 
®c22 ^c33 difficult. Thus, try a three section lead 


compensator . 
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4) By spirule (or by inspection), three zeros close to the origin 
and two poles far away from the origin can stabilize Figure 8.2 
to Figure 8.5 to a great extent. For one specific choice. 


ki(s + n Cs + 0.S)2 

s(s+4) (s+10) 


where is a free parameter. 


the results are shown in Figure 8.6 to Figure 8,9, Note that 
the loci for kj<0 are not shown, since k^<0 is not desirable 
for the same reason as in 1 (i) . 


5) By inspection of these root i^oci and by considering the effect 

of kj on the pole-zero pattern of Fj and ^3/^2 Csee C7.36) and 
(7.39)), kj^ = 1000 is chosen for the same reason as in the 
previous design of Section 7.3, For this value of k^, the 
root loci for F^ and F3/P2 drawn. Then, g ^.22 

designed according to these two root loci and their relationship 
with the poles and zeros of G (see (741)). By trial-and-error 
(or by inspection) , gc 22 ~^^ (s+O . 5) ^ free parameter) 

was found to be good. For this 8 ^. 22 * root loci for (7.16) 

and (7,17) are shown in Figure 8.10 and Figure 8 , 11 . Loci for 
k 2>0 are not shown since for k 2 > 0 , one branch in each plot 
extends along positive real axis to + hence undesirable. 

6 ) By inspection of Figure 8.10 and Figure 8.11 and by (7.41), 

IC 2 « “400 is a good value. 










115 


7) Using k 2 = -400, the root locus for (7.13) with g ^.32 

s ^3 (s-«-0.3) (s + 0,4) drawn (see (7,41)), and g ^33 can be 

s (s+10) 

designed accordingly. The result is shown in Figure 8.12, from 
which the best value for can be seen by inspection to be 
k3=-20. 

Thus, a design with 

gcii(s) = I000(s+l)(s+0.5)^/s(s+4)(s+10) 

" -400(s+0.5)/s (8.2) 

gc 33 (s) = -20(s+0.3) (s+0.4)/s(sfl0) 

is completed. 

The schematic diagram for the designed system is shown in 
Figure 8.13* 

For step inputs of magnitudes 126.7 ft, /sec,, -0.25 radian, -0.5 
radian in tj^,r 2 » and r^, the simulation results for y 2 , and 

are shown in Figure 8,14, Figure 8.15 and Figure 8.16, 
respectively. It is seen that both stability and steady-state 
decoupling have been achieved. Furthermore, due to the introduced 
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pole at the origin in each of ^c33' steady- 

state error in each of the outputs is zero, which is the most 
desirable situation. Thus, by, decoupling the steady states, 
steady-state accuracy has also been achieved as a byproduct. 



FIGURE 8.13 






i.OQ a. 00 !?■ 00 16 . 00 20 . CO 24.03 

FIGURE 8.14 

TIME RESPONSE OF OUTPUT IN FIGURE 8,13 WITH 


rj=126.7'u(t) 
r2=“0,25*u(t) 
rj=-0.5’u(t) 
u(t) : UNIT STEP 



FIGURE 8 . 15 


TIME RESPONSE OF OUTPUT IN FIGURE 8.13 WITH 

r^=126.7-u(t) 

r2=-0.25’u(t) 

, 5 ' u (t ) 


uCt) : UNIT STEP 



ii 


FIGURE 8.16 

TIME RESPONSE OF OUTPUT IN FIGURE 8 

r^sl26.7-u(t) 
r 2 = -0. 25 *u(t) 
r3*-0.5‘u(t) 
uCt) : UNIT STEP 


13 WITH 
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9, CONCLUSIONS 

9.1 COMPARISON OF THE RESULT IN PART I TO THAT OF STATE 
VARIABLE FEEDBACK APPROACH 

In part I, a constructive criterion for decoupling the steady states 
of a linear time-invariant multivariable system was developed. 
Transfer function matrix representation, unity feedbacks and 
cascade compensation were used as shown in Figure 2,1, 

Another approach using linear state variable feedbacks was 
investigated by Wolovich (2) . The result in terms of transfer 
function matrix representation is given in Chapter 1, and repeated 
here as follows: 

A system characterized by an (nxm) proper rational transfer 
function matrix, G^Cs), having no poles at the origin 
(s=0) can be steady-state decoupled (via linear state 
variable feedback or perhaps some other less ambitious 
scheme) if and only if 

p(Gp(0)) -n (9.1) 

where p(Gp(0)) denotes the rank of the matrix G^Cs) when s 
approaches zero . 
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Several advantages of the result given in Part I over that described 
above are in order: 

1. Cascade compensators and output feedbacks are much easier to 
implement than measuring the states, 

2. The rank condition (9.1) is not necessary, e.g,, given the 
2 x 2 plant (n=m = 2 ) ; 

1 

s + 2 
1 

s+ 4 

By (3.19), (3.10) and (3.18), steady states can be decoupled 

by introducing two and one pure integrators in g^. and 8^22 
in Figure 3.1, respectively. However, by (9.1), this cannot 
be done through linear state variable feedbacks. 

3. Poles at the origin in the given plant are allowed. Actually, 
such poles are very helpful for steady- state decoupling as 
was shown in Figure 6.6 for the plant (6.24). 

9.2 POSSIBLE GENERALIZATION OF THE SYSTEM CONFIGURATION 

The discussion of this thesis has been restricted to the configur- 
ation given in Figure 2.1. For the general feedback configuration 
in which the unity feedbacks are replaced by a transfer function 
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matrix G^Cs) (nxn) , the simple relation (2.5) is no longer valid. 
However, if G^(s) is diagonal and nonsingular (in the field of 
rational functions of s), which are of practical importance, similar 
result to (2.5) can still be obtained as follows: 

H = (I*GpG^Gf)-lGpG^ 

= (I.GpGcG^)-^GpG^GfGf-l 

= (I-(I.GpG^G^r^) 

Since is diagonal, is also diagonal. Therefore, simple 

expressions for the off-diagonal elements of H(s) can still be 
obtained easily. Then, with a slight modification, the results in 
Part I can still be applied. 

9.3 STABILITY AND DESIGN 

In Parts II and III, stability of a linear time-invariant multi- 
variable system was considered. A design technique, using an 
extended root locus method was also developed and applied 
successfully to 2x2 and 3x3 cases. The major achievement is the 
revelation of the simple connection between single-loop and 
multivariable cases. Such connections made the application of 
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singl e- loop design methods to multivariable systems possible, as 
was seen through the design examples in Chapter 6, Chapter 7 and 
Chapter 8. Some other advantages of the design techniques are: 

1. Consideration of integrity problems is possible in the 

process of the design by forming pertinent root loci. 

This means that the system can be designed such that possible 
failure of any loop (or combination of loops) do not cause 
the system to be unstable Ce.g., see (10)). 

The problem of input output permutation, like ’’which output 
should be fed back to a particular input?” can be solved to 
some extentbyinspection, 

3, More insight to the problem is achieved through the root 
locus approach . 

4. The problem of meeting system specifications can be done in 
the same manner as in any single-loop design method. 

A few disadvantages, however, do exist. For example, the successive 
dependence of each root locus on the previous ones causes more design 
difficulty as the number of inputs and outputs of the plant increases. 
Also, like the single-loop frequency domain methods, trial -and -error 
is inherent in this technique. 
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However, with the help of computers , these problems can be minimized 
and the design can be done within a reasonable amount of time. 
Besides, with some experience in handling the root locus, the effort 
can be further reduced. 



12S 


REFERENCES 


(1) P, L, Falb and W, A. Wolovich, ’’Decoupling In The Design 
and Synthesis of Multivariable Control Systems” IEEE 
Trans. Auto. Control Vol. AC-12(6), December, 1967, 

pp . 65 1 - 659 . 

(2) W. A. Wolovich, ’’Static Decoupling” Joint Auto. Control 
Conference, 1973. 

(3) E, G. Gilbert, ’’Controllability and Observability in 
Multivariable Control Systems” SIAM J, Control Ser, A 
Vol . 2, pp. 128-151, 1963. 

(4) C. T. Chen, "Representations of Linear Time- Invar iant 
Composite Systems” IEEE Trans. Auto. Control, Vol, AC-13(3), 
June, 1968. 

(5) H. H. Rosenbrock, "Design of Multivariable Control Systems 
Using the Inverse Nyquist Array" Proc . lEE, Vol. 116, 

pp. 1929-1936, 1969. 

(6) F. R. Gantmacher, "The Theory of Matrices" Vol. 1, pp. 8-12, 
Chelsea Publishing Co., 1960, 

(7) C. T. Chen, "Stability of Linear Multivariable Feedback 
Systems” Proc. IEEE, Vol, 56, No. 5, May, 1968, 

(8) C. H, Hsu and C. T. Chen, "A Proof of The Stability of 
Multivariable Feedback Systems" Proc. IEEE, Nov., 1968. 

(9) C. T. Chen and C, H. Hsu, "Design of Dynamic Compensators 
for Multivariable Systems" Proc. 1971 JACC, pp , 893-900, 

(10) J. J. Belletrutti and A. G. J. MacFarlane, "Characteristic 
Loci Techniques in Multivariable-Control-System Design" 

Proc. lEE, Vol. 118, pp . 1291-1297, 1971. 

(11) A. G. J. MacFarlane, "Return-Difference and Return-Ratio 
Matrices and Their Use in Analysis and Design of 
Multivariable Feedback Control Systems” Proc. lEE, Vol. 117, 
pp. 2037-2049, 1970. 

(12) B. D. 0. Anderson, "Linear Multivariable Control Systems - 
A Survey” Proc. IFAC 1972. 



126 


(13) 

(14) 

(15) 


C, T. Chen, "Introduction to Linear System Theory" Holt, 
Rinehart and Winston, New York, 1970. 

G. J. Thaler, "Design of Feedback Systems" Dowden, 
Hutchinson § Ross, Inc., 1973. 

D. Hazony, "Elements of Network Synthesis" Reinhold, 

New York, 1963, 



127 


APPENDIX A: FORMULAE (4.1) AND (4.2) 


The proof of (4.1) can be done by applying Lemma A.l in Appendix 
A of (15) directly to the two nxn square matrices I and G. For 
(4.2) however, direct application of this lemman has a little 
difficulty. 

For better analytic ity , an independent proof using mathematical 
induction has been developed. An outline of the proof is given 
below: 

1. (4.1), (4.2) are satisfied for n=2 and n=3 by direct 
expansion. 

2. Suppose (4.1) is true for n^N, (4,2) is true for n=N+l, 
then (4.1) is true for n=N+l. 

3. Suppose (4.1) is true for both n=N-l, and n=N, then (4.2) 
is true for n-N+1. 

Thus, starting from N*2, it can be induced that (4.1) and (4.2) 
are true for any positive integer. 

Details of the proof in 2 and 3 above are omitted. 
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APPENDIX B; CHARACTERISTIC EQUATION 


It was pointed out in Section 6.1 that the stability of a multi- 
variable system as shown in Figure 2.1 is determined by the zeros 
of both (6.2) and (6.3). It will be shown in this appendix that 
if pole-zero cancellations are done deliberately, and if the 
compensator poles are chosen carefully, zeros of (6.2) alone can 
determine the stability. The constraints under which this is true 
are very practical and can be fulfilled in a systematic manner* 
Thus, the mathematical possibilities in which zeros of (6.3) must 
be considered are bypassed. 


By the definition of characteristic polynomials given in Section 
6.1, and dp(s) can be expressed analytically as: 


A.(s)*LCD{q,| 


/il.--*,i 


lsjj< 


min(n,m), liij^<***<ij^^m, 
<jj^<n) (B.l) 


(s)=LCD (G 



je,*l , * • • ,min (n, m) , 
Isj 1< ‘ ‘ ’ <3 ^ 


lsij<’ • *<i^<n, 

(B.2) 


Where, as in Chapter 4, 
% 




denote the 
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£th-order minors of G and G formed from rows ii»** 

c p 

columns 3i»'’’»3£ of each matrix, respectively, 

LCD {— •} denotes the Least Common Denominator of all the rational 
functions described in the brackets and all the minors are 
assumed to be in irreducible rational forms. 




denote the numerators and denominators of the irreducible 




respectively . 


Then, 


A^.(s)=LCM{D^. 



I £ = 1 , * • • ,min (n,m) I5 
1< j * • •< j^<n} (B.3) 



(s)=LCM{D 



il, ‘ 

3 l»‘ 



|£ = 1,’ • *,min(n,m) l<ij^<»**<i^<n, 

lsjj<- • (B.4) 


Where LCM {••*} denotes the Least Common Multiplier of all the 
polynomials described in the brackets. 
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These analytical expressions for ^^,Cs) and ^p(c) will be useful 
later in this appendix. 


For single-loop systems. 


Gp(s) = [gp(s)]j^^ 

G^(s) = [g^(s)]i^l 


Let N (s) , N (s) , D (s) , 0 (s) denote the numerators and the 

U v« P w 

denominators of g„(s) and g^(s) (both in irreducible rational forms) 

P ^ 

respectively, we have for (6.2) 

Ni(s) N (s) N^.(s) 

oTTiT " ^ " Dp(s) ■ D^(s) 

(Note that Nj^(s)/Dj(s) is in irreducible form.) 

Also, by definition of the characteristic polynomial, 


fip(s) = Dp(s) 
A^CsD = D^(S) 



Hence, (6,3) becomes 


N(s} 


D^Cs) DpCs) 
DiCs) 


The right hand side of (B.5) can be written as 


(B.6) 


Dp(s) D (s) + N (s) N^,Cs) 

— Dp(s) t)V(H ^ 

Let A 

Nq(s) a DpCs) D^(s) + Np(s) Nj,(s) 

which are the denominator and numerator of (6,2) before cancellation 
(if any). Also, let C(s) denote the greatest common factor between 
D^(s) and N^Cs) CC(s)= 1, if no common factor exists), then 

D„(s) = C(0 • DjCs) 

n^(s) = C(s) ■ NjCs) 


By (B.6) and (B.7), we have; 
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(B,8) tells us that the zeros of Nj(s), together with those pf N(s) 
are the zeros of N^(s) alone. 

Therefore, the following conclusion, which is well-known in single- 
loop theory, can be made. 

The stability for systems shown in Figure B.l is determined by the 
roots of 

1 + gp(s)g^(s) = 0 (B.9) 

if and only if no pole-zero cancellation is allowed in (B.9), even 
if a common factor exists. 


■^) 

gcCs) 


gpCs) 











FIGURE B.l 


As an example, consider the system shown in Figure B.l with 


gp(s) 


= Izl 

s + 1 




1 

s- 1 
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By (6.1) and (6.3), 


Nj(s) 


S + 2 

Dj(s) 


s + 1 


A 


(s-lHs+l) 

N(s) 

s 

— in 


A , 

Since the zero of N(s) is in the right-half plane, the system is ? 
unstable. 


By (B.9), 




0 


s + 2 

If the common factor (s-1) is cancelled, we have j^-0. Therefore, 

♦ 

only thes=-2 pole is retained and erroneous conclusion that the 
system is stable is reached. 

However, if Cs-1) is not cancelled, we have Viy = ^ Q > Hence, 

both the zero of N^(s) and N(s) are retained. 

Therefore, for single-loop systems, if no cancellation is allowed, 

A ' ' 

(6,2) alone gives all the zeros of Nj(s) and N(s), hence determines 


the stability. 
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The following question then arises naturally: ’’Can we use (6.2) 

alone in determining the stability of a multivariable system by the 
same requirement that no cancellation is allowed?” The answer is 
”no” as was shown by Chen (7) through the following 2x2 example; 


Consider 


G (s) 
P 


“S^+S+1 
(s+T)'^ (s-l) 

1 

Ts+i) (s-i) 


1 

s-l 

1 

s-l 


(B.IO) 


G^(s) = I 


(B.ll) 


where I denotes the 2x2 unity matrix. 


Then, 


det(I + G G ) = s 1 


= cs-n 

(s-l)2(s + l) 


(B.12) 


There exists a right-half-plane zero in (B.I2) at s*!. 


However, by (6.2) and (6.3) 
Nj(s) * 1 


« s + 1 


(B.13) 


N(s) 
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(Note^A^Cs) = 1, Ap(s) = (s+l)(s-l) and D^(s) « s-1) the system is 
clearly stable. Thus, leaving all the common factors uncancdled 
does not work , 

However, if the general formula (4,5) is used and let all the minors 
of G be in irreducible forms, we have 

P 

det(I+G^G^) = det(I+G^) 
pc p 

®pll ^P22 p 

_i^ -s^+s+1 . 1 ^ -s 

' ^ TsTryxFTT s^ cs+i)‘(s-iT 

s + l 

' (s + nTs-TT . (B.14) 

It is seen that if the common factor (s+l) is not cancelled, the 
zero of (B.14) is exactly the same as those given by Nj^(s) and 

A ’ 

N(s) as shown in (B.13). 

Note that the misleading factor (s-1) in (B.12) does not appear in 
(B.14). The reason is that we started out with irreducible minors 
and in forming the 2nd order irreducible minor detGp, the (s-1) 
factor was cancelled. 
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Therefore, by selecting cancellations, stability of a multivariable 
system can still be determined by (6.2) alone. All common factors 
in the minors must be cancelled to get the irreducible forms, while 
the others are not allowed. This strict rule together with the 
application of formula (4.5) make the whole procedure completely 
systematic, no confusion will arise. 


As another example, consider 



G (s) = I 
c 

which is also an example in (7) . 

By direct manipulation, it can be found that det (I+G G ) ® -J. 
However, by selecting cancellations as described above, we have 



(B.15) 
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Once again^ if the common factor (s + 1) (s-1) is not cancelled, the 
two zeros of (B.15) are exactly what one would obtain as zeros for 
Nj(s) and N(s) by (6.2) and (6,3). Thus, these two uncancelled zeros 
of det(I+GpG^) determine the stability of the system. Since one of 
them is in the right-half -plane , the feedback system is unstable. 


These two examples suggest that (6.2) alone can be used as 
characteristic equation for a multivariable system if we select 
cancellations as described above. But, is this true in general? 

To answer this question, consider the general expression (4.5) 

for feedback system as shown in Figure 2.1. For simplicity, consider 

2x2 and first. 

By (4.'5), 


det(I^GpG^) = l*Gp{J)GjJ)*Gp(l)G^(2).Gp(2)G^(l).Gp(2)G^(2) 

CM, 2 


*e 1 + 


"jj) Np(^) 


v\:P 


(B.16) 


Where the notations for numerator and denominator of each irreducible 
minor used in (B.3) and (B.4) are employed, e.g., Np(J), Dp(J) 
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denote the numerator and denominator of the first order minor 
Gp(J) * etc. Note that no common factor exists between 

Np(l) and Dp(J). N^(J) and D^(l), etc. 

Let N^(s), D^(s) be the numerator and denominator of (B.16), after 
collecting all the terms at the right-hand-side under the 
restriction that no pole-zero cancellation is allowed (even if a 
common factor exists). Then, 

Dj,(s) = LCM{Dp(i)D^(J),Dp(l)Dc(f),Dp(J)Dg(^), 

Also, by (B.3) and (B.4), 

ip(s)=LCM{Dp(J).Dp(l).Dp(2).Dp(2),Dp(J.2)} (B.19) 

If no common factor exists between ^*(s) and ^ (s) , which can be 

c p 

realized by not using any plant pole^ as a pole in the compensator, 
^G(s) has a pole at s=X, whenever at least one element of G(X) is ~. 
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we have 


CB.20) 


This can be proved by the following arguments: 


1 . 


If no common factor exists between A-Cs) and A (s), then no 

p 

common factor can exist between any element in (B,18) and any 

of those in (B.19) Otherwise, a common factor will exist 

between A^(s) and A_Cs). 

^ P 


2. Any factor of must be a factor of either A^(s) or 

Ap(s) and with the same multiplicity. The reason is that any 

factor of Dq(s) must exist in at least one of the five 

elements in (B.17). And the multiplicity of this factor 

must be the same as that of the element that has the maximum 

multiplicity of the same factor. By 1, this factor can be 

either in the D *s or in the D *s of (B.17), but npt both. 

P 

Therefore, by (B.19) or (B.18), the factor with the same 

multiplicity must appear in either A„(s) or A_(s). 

P c 

3. Any factors of A^(s) must also be a factor of D^Cs) 

with the same mult iplicity. This can be seen by similar 
arguments as that in 2 above, but starting from (B.18) and 
(B.19) instead of (B.17). 
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Now, similar to what was done in the 1x1 case, let C(s) be the 
common factor between N (s) and D^jCs). Then, 

I>o(s) ' C(s)-Di(s) 

N^(s) = C(s)-NiCs) 

By (6.3) and (B.20), 

N(s) = ^cCs) 

J>lls) 

« C(s) 

Therefore, the zeros of N^(s) are exactly those of N^Cs) together 

' A , 

with those of N(s), Thus, stability of a 2x2 multivariable system 
can be considered by (6.2) alone; If (1) irreducible minors are 
used in (4.5), (2) all the other cancellations are not allowed, 

(3) no plant pole is used in the compensator, and (4) g . . (s) ^ 0 

Clj 

for alli=l,'’*,mjal,***,n. 


The condition gpi*(s) 0 in (4) was added because any zero c 
will cause one corresponding element of (B.17) to be missing which 
will impair the equality of (B.20) if the associated Dp term happens 
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to be the only one among all the Dp’s that contain the highest 
multiplicity of any factor. This is shown by the following example 


r 1 


Let 


Gp(s) 


iTT 

s-5 

(s + 5)(iVl) 


1 

7F3T 


s + 2 


GJs) = 


gcii(s) 

8c2iCs5 


*c22 


Then 


Kpll^s) 


s + 1 






S-5 

TTnynny 


g_ _(s) 
P22 


s + 2 


detG^(s) « 1 

P “(s+i) (s + 2) (s + 3) 


(B.21) 
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By (B . 19) , we have 

Ap(s) = LCM{(s + l), (s-5), (s + 3)(s + l), Cs + 2), (s + 1) (s + 2) (s + 3) } 

= (s+1) (s+2) (s+3) (s-5) 

Note that the highest multiplicity for the factor Cs-5) comes from 

the second element in the bracket, which is D (i) , denominator of 

P ^ 

^c21^^^ = 0 is used (e.g., diagonal G^,) , 
the third term in (B.16) is zero, hence the element D (i) D (?) will 

P C i 

not appear in (B.17). Therefore, the factor (s-5) will not appear 
in D^(s). This makes (B.20) to be not true. The (s-5) factor will 

A 

show up in N(s) of (6,3) which makes the system always unstable as 
long as Therefore, it is impossible to stabilize the 

system with diagonal compensator matrix G^(s). 

Therefore, whenever diagonal G^(s) is employed, care must be taken 
to see if situation like this happens. This can be checked very 
easily by forming all the pertinent irreducible minors of Gp(s) and 
then check to see if the multiplicity of each plant pole is retained 
in these minors. If yes. zero element is allowed in GJs) and 
(6.2) alone can determine the stability. If not, (6.3) must also 
be considered as was shown in the above example. Incidentally, this 
is a trivial case in single-loop systems. Since there is only one 
element in G^ (s) , which obviously cannot be identically zero. 
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The above analysis for 2x2 Gp(s) can be carried over in exactly the 
same manner to the nxm case. 

Therefore, for systems as shown in Figure 2.1, if 

1. The general formula (4.5) for det(I+G_G ) is used. 

Pc 

2. All the common factors in the minors of Gp are cancelled to 
get irreducible forms. 

3. All other cancellations are not allowed. 

4. Plant poles are not used as poles of G^ , 

$. Appearance of zero elements in G^, is carefully checked, 

then the zeros of detCl + GpG^,) alone determine the stability of the 
system. If all of them are in the open left-half-plane, the system 
is stable, otherwise it is not . 

Therefore , as in the single- loop theory, 

det(I+GpGg)=0 (B.22) 

is referred to as the characteristic equation for multivariable 
systems and stability can be considered through this equation. 
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APPENDIX C; ROOT LOCUS GAIN 


When root locus approach is used, it is convenient to express 
transfer functions in terms of root locus forms ( 14 ) as shown 
below: 


G(s) 


kJCs+z j) 

|jk 

s n(s+pj^) 


(c.i) 


where j, k are positive integers and s=-Zj, s=-pj^ are the zeros 
and poles of G(s) respectively. 

The gain constant k in (C.I) is referred to as the ^'root locus 
gain'* of the function G(s) (see 14 ). 

In step 6 of Section 6 . 4 , some algebraic manipulation was performed 

to find the root locus gain k__ of the function G (s) . In most 

eq 

cases, this step can be bypassed as shown below: 


Let kj^, k2 be the root locus gains of the two transfer functions 
Gj(s) and 62(5) respectively. Let CGg(s), G^(s) be the sum and 
the ratio (G2/Gj^) of Gj^(s) and 02(5), in root locus forms and k , 
k^ denote the corresponding root locus gains. Also, let the order 
of the numerator and denominator of any rational function H(s) be 
denoted by ON(H) and OD(H) respectively . 



Then, we have 


THEOREM C. 1 

The root locus gain of G (s) =G^ (s) +0^ (s) is 

(i) kg = kj if and only if ON (G j ) +0D (G2) >0N (03) +0D (G^) 

(ii) k^ = k2 if and only if ON (G^) +0D (G2) <0N(G2) +0D (G 

(iii) k^ = ^ 1+^2 ON(Gj)+OD(G 2)=ON(G2)+OD(G^) 

and ^ -k2 

THEOREM C.2 

The root locus gain k^ of G(s) = G2(s)/G^(s) is 

Proofs for both Theorem C.l and C.2 are straightforward, hence, 
omitted. 

Both of these theorems are very simple in nature, however, they are 
very useful tools in evaluating the root locus gains, as illustrated 
below for the determination of k in Section 6.4 (see (6.41)). 
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By (6.8), 

= l*Bpll8cll 

and G. = 2l 

eq Gi 

For G2 » 

OH(gp22^*OD(detGp-g^ll-) = 0*4 = 4 
ON(detGp-gpll)+OD(gp22) -2+1=3 

Since 4>3, we have by Theorem C.l (i) 



where . k denote the root locus gains of G2^®^ 

2 ®p22 ^ 

respectively. 


Similarly, for Gj^ 

0N(l)+0DCgpj^gj.ll) = 0 + 2 = 2 

0N(gp^jg^ll)+0D(l) = 1+0 = 1 
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Again, since 2>1, by Theorem C.l (i), we have 



Then, by Theorem C.2, 



which agrees with what was obtained in Section 6.4 through algebraic 
manipulation. 

* 

Remark: Whenever kj+k2*0 in case (iii) of Theorem C.l, no conclusion 

can be obtained through the theorem , since no obvious analytic 
expression exists for the coefficient of the second higher order 
terms. However, direct algebraic manipulation can always be used 
in such cases. 



